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SUMMARY 


This  report  describes  the  work  performed  under  Contract 
AF  04 (611) -9073  during  the  second  quarter  of  the  contract  period, 
beginning  September  1963.  The  objectives  of  the  study  are: 

(a)  To  develop  design  tools,  in  the  form  of  digital  computer 
programs,  for  the  nozzle  designers  to  use  for  calculating  the 
behavior  of  materials  exposed  to  a  solid-propellant  rocket-exhaust 
environment . 

(b)  To  marshall  existing  theoretical  methods,  and  to  extend 
these  or  develop  new  methods  where  necessary,  to  gain  a  better 
understanding  of  the  fundamental  mechanisms  associated  with  exposed 
materials  behavior  in  a  rocket  nozzle  environment. 

The  work  during  this  quarter  has  been  directed  toward  the 
development  of  a  digital  computer  program  applicable  to  materials 
such  as  tungsten  and  graphite,  that,  in  general,  do  not  erode 
appreciably  or  that  erode  with  one  moving  boundary.  This  has 
required  effort  in  the  areas  of  convective  heat  transfer,  transient 
conduction,  equilibrium  chemistry,  reaction  kinetics,  and  mass 
transport,  as  applied  to  the  solid-propellant  rocket  problem,  and 
this  work  is  discussed  in  detail.  Further,  a  plasma-generator 
rocket  simulation-test  program  has  been  formulated  to  complement 
the  theoretical  studies,  and  the  test  program  and  some  firing 
results  are  described. 
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A  STUDY  OF  SOLID-PROPELLANT  ROCKET 
MOTOR  EXPOSED  MATERIALS  BEHAVIOR 

I.  INTRODUCTION 

During  the  past  several  years,  significant  theoretical  tech¬ 
niques  have  been  developed  by  various  investigators  for  character¬ 
izing  many  of  the  complex  individual  phenomena  occurring  in  the 
ablation  process  of  solid-propellant  rocket  wall  materials.  These 
techniques  have  been  used  to  gain  a  better  understanding  of  abla¬ 
tion  phenomena  under  simplified  conditions  and  to  guide,  to  a 
certain  extent,  the  design  of  aft-closures  and  nozzles.  The  proven 
utilitiy  of  the  theoretical  approach  in  this  area  has  made  meaning¬ 
ful  the  possibility  of  developing  an  integrated,  more  general,  ana¬ 
lytical  technique  for  the  design  of  solid-propellant  rocket  wall 
materials  and  configurations.  The  need  for  such  a  technique  is 
becoming  more  urgent  as  motors  increase  in  size,  and  thus  render 
the  traditional  cut-and-try  approach  more  costly  and  time-consuming 

Recognizing  this  need,  the  objectives  of  this  study  are  two¬ 
fold: 

(1)  To  develop  design  tools  for  the  nozzle  designer  to  use 
for  calculating  the  oehavior  of  materials!  exposed  to  a  solid- 
propellant  rocket-exhaust  environment.  These  design  tools  will  be 
in  the  form  of  digital  computer  programs. 

(2)  To  marshall  existing  theoretical  methods,  and  to  extend 
these  or  develop  new  methods  where  necessary,  to  gain  a  better 
understanding  of  the  fundamental  mechanisms  associated  with  exposed 
materials  behavior  in  a  rocket  nozzle  environment.  This  second 
objective  is,  of  course,  coupled  with  the  first  in  that  an  appre¬ 
ciation  of  the  fundamental  ablation  mechanisms  is  required  to  calu- 
late  their  effects. 

Both  theoretical  and  experimental  studies  are  being  conducted 
to  accomplish  the  objectives  of  the  program.  The  relations  charac¬ 
terizing  the  ablation  process  are  being  formulated  largely  from 
theoretical  considerations,  but  also  draw  on  experimental  work  in 


those  areas  not  amenable  to  theoretical  treatment.  Existing  test 
data  are  being  used  where  avilable,  and  small-scale  laboratory 
tests  are  being  conducted  as  a  part  of  this  program  to  provide 
information  in  several  areas  where  it  is  required.  The  flow  chart 
of  Figure  1  indicates  how  the  various  areas  of  investigations  are 
interrelated . 

As  indicated  in  Figure  1,  the  theoretical  studies  are  divided 
into  two  major  categories,  the  first  having  to  do  with  materials 
that  erode  with  one  "moving  boundary"  such  as  graphite  (nonablating 
refractories  are  a  special  case  of  this  category) ,  and  the,  second 
with  materials  such  as  the  reinforced  plastics,  or  composites, 
whose  ablation  can  be  characterized  by  two  moving  boundaries.  The 
experimental  studies  are  being  carried  out  in  the  Vidya  1-megawatt 
plasma-generator  facility;  they  consist  of  graphite  and  composite- 
nozzle  ablation  tests  and  metallic-oxide  particle-impaction  tests. 

( 

Durirg  the  second  quarter  of  the  program,  the  theoretical 
effort  continued  in  the  following  areas:  convective  heat  transfer, 
thermal  behavior  of  a  nonmelting,  noncharring  nozzle  insert,  nozzle 
wall  chemical  erosion,  and  reaction  kinetics.  The  accomplishments 
here  are  reported  in  Section  2.  The  experimental  effort  during  the 
quarter  consisted  of  facility  modification  and  checkout  for  the 
required  operating  conditions,  test  section  fabrication,  calibration 
firings,  and  an  ATJ  graphite  nozzle  firing.  This  work  is  discussed 
in  Section  3. 

2.  THEORETICAL  STUDIES 

2.1  Introduction 

For  the  purpose  of  classification,  all  engineering  effort  on 
the  program  that  does  not  directly  involve  an  experimental  test 
program  is  called  "theoretical."  This  is  actually  a  misnomer,  since 
the  effort  is  based  largely  on  proven  physical  and  mathematical 


relations,  and  is  of  a  very  practical  and  applied  nature.  That 

is,  the  effort  is  directed  toward  the  construction  of  a  calculation 
tool  for  characterizing  the  phenomena  occurring  at  an  ablating 
nozzle  wall.  Because  many  physical  relations  are  needed  for  this 
characterization,  and  because  the  interactions  between  them  are 
quite  complex,  an  analysis  flow  diagram  was  prepared  (Fig.  2)  to 
help  clarify  the  role  of  each  area  of  study. 

Figure  2  is  itself  necessarily  complex,  but  it  is  hoped  that 

it.  will  be  helpful  in  placing  each  phase  of  the  effort  discussed 
herein  in  proper  perspective.  The  row  of  boxes  at  the  top  of  the 
figure  represents  the  various  inputs  required  to  perform  nozzle  wall 
erosion  calculations;  except  for  the  erosion  effects  of  particle 
impact,  none  of  these  specific  data  are  being  sought  experimentally 
under  the  present  contract.  The  program  emphasis,  rather,  is  on 
the  "calculation"  area  of  the  figure,  to  provide  a  technique  for 
determining  nozzle  erosion  rates  and  wall  temperature  from  the 
input  data. 

The  work  reported  in  this  section  is  in  the  areas  of  transfer 
coefficients,  surface  chemical  reactions,  and  internal  material 
behavior.  As  the  reader  reviews  the  work,  it  is  recommended  that 
he  refer  to  Figure  2  occasionally  to  review  the  role  of  each  study. 

2.2  Boundary-Layer  Transfer  Coefficients 

The  boundary- layer  transfer  coefficients  required  for  the  cal¬ 
culation  of  nozzle  wall  recession  rate  are  those  of  heat,  C„ , 

ii 

mass,  CM,  and  momentum  C^,.  These  three  coefficients  are  actually 
interrelated,  and  their  exact  determination  would  follow  from  a 
simultaneous  solution  of  the  boundary- layer  energy r  diffusion,  and 
momentum  equations.  It  has  become  common  practice,  however,  to 
employ  in  their  determination  rather  simple  approximate  relations. 
Along  this  line,  the  approach  used  herein  is  to  calculate  C„  inde- 

n 

pendently  of  CM  and  CF ,  then  to  estimate  from  it  using  the 


Chilton-Colburn  relation  (although  the  validity  of  this  requires 

examination  for  the  multi-component  boundary- layer  situation) . 

The  momentum  transport  coefficient,  Cp ,  is  not  required,  except 

when  dealing  with  materials  susceptible  to  mechanical  erosion  due 

to  wall  shear;  in  these  cases,  it  should  be  adequate  to  determine 

C_  from  C„  with  Reynolds'  analogy.  Hence,  the  major  effort  to 
r  ri 

date  has  been  in  the  selection  of  a  technique  to  calculate  C„. 

ri 

2.2.1  Nonablating  nozzle  wall  convective  heat  transfer 

Calculation  of  the  convective  heat-transfer  coefficient  for 
a  nonablating  nozzle  wall  is  an  important  first  step  in  the  deter¬ 
mination  of  nozzle  ablation  characteristics,  since  this  value, 
after  correction  for  the  blowing  reduction  (a  small  correction  in 
this  case) ,  is  an  integral  part  of  the  wall  ablation  relations  as 
shown  in  the  First  Quarterly  Progress  Report,  Reference  1.  A 
survey  discussion  of  the  several  convective  computational  methods 
was  presented  in  that  report,  and  it  was  pointed  out  that  the  ARM 
method  (after  Ambrok ,  Rubesin,  and  Mayer)  appeared  to  offer  the 
advantages  of  reasonable  accuracy  and  simplicity  of  utilization. 
During  the  past  quarter  the  ARM  method  was  compared  with  experi¬ 
mental  data  obtained  at  the  Jet  Propulsion  Laboratory,  and  the 
results  were  indeed  encouraging.  Although  more  checks  will  have 
to  be  made  (particularly  with  rocket  data)  before  the  method  can 
be  given  an  unqualified  recommendation,  it  appears  at  this  junc¬ 
ture  that  it  is  superior  to  others  in  common  use.  Hence,  its 
development  is  given  in  detail  herein,  and  a  comparison  of  its 
results  with  experimental  data  is  presented. 

2. 2. 1.1  The  Ambrok-Rubesin-Mayer  Method 

The  essence  of  the  ARM  method  is  the  assumption  that  the 
flat-plate  relationship  between  the  local  boundary- layer  energy 
thickness,  0,  and  the  Stanton  number,  C„ ,  is  valid  also  for 

rl 


flows  with  a  streamwise  pressure  gradient.  This  assumption  permits 

solution  of  the  pressure  gradient  boundary- layer  energy  integral 

equation  for  0,  and  then  employed  again,  permits  calculation  of 

the  pressure  gradient  Stanton  number,  C„,  from  0.  This  is 

£1 

explained  in  greater  detail  later  in  the  section.  For  comparative 

purposes,  the  alternate  approach  in  most  common  usage  today  should 

also  be  Mentioned.  It  is  the  use  of  the  boundary- layer  momentum 

integral  equation,  wherein  the  flat-plate  relation  between  momentum 

thickness,  0,  and  friction  factor,  Cp ,  is  assumed,  along  with 

two  additional  assumptions:  (1)  velocity  profile,  to  permit  calcu- 

* 

lation  of  the  shape  factor  5/0  (appearing  in  the  momentum  equa¬ 
tion)  ,  and  (2)  Reynolds  analogy,  to  permit  calculation  of  C„  from 

n 

the  resulting  Cp. 

The  historical  development  of  the  ARM  method  is  of  some 
interest,  inasmuch  as  it  provides  some  further  insight  to  its 
validity.  To  the  best  of  the  writer's  knowledge,  the  basis  of  the 
idea  developed  first  in  this  country  in  the  early  1950's  when 
Rubesin  (Ref.  2)  at  the  Ames  Research  Center  of  the  NACA  discovered 
that  convective  heat-transfer  measurements  could  be  correlated  quite 
well  with  the  local  boundary- layer  energy  thickness.  Later,  at  Vidya, 
Rubesin  used  this  knowledge  to  derive  the  present  method,  and  first 
published  it  in  Reference  3.  The  application  of  interest  at  that 
time  was  aerodynamic  heating,  and  it  has  been  in  use  at  Vidya  for 
3  years  for  that  purpose.  Its  results  have  been  compared  favorably 
with  re-entry  body  and  X-15  convective  heating. 

Independently,  Ambrok  suggested  the  essence  of  the  same 
technique  (Ref.  4) ,  but  his  paper  was  not  noted  by  the  workers  at 
Vidya  until  Mayer  (Ref.  5)  presented  an  application  of  it  to  the 
rocket-nozzle  situation.  Since  it  is  felt  that  the  contributions 
of  all  three  investigators  are  important,  the  appellation  "ARM"  is 
used  for  the  method  herein. 

The  derivation  of  the  ARM  method  follows  from  the  integral 
form  of  the  boundary- layer  energy  equation  (with  no  blowing) : 
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where 


and 
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p'u  (H  -  H  ) 
r  e  r  w 
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0  = 


£U 


dy 


(3) 


Note  that  this  latter  definition  differs  from  the  conventional 
energy  thickness  definition  in  which  is  used  in  place  of  HrJ 

this  modification  is  required  here  to  prevent  an  unnatural  singu¬ 
larity  in  the  resulting  equations  in  the  case  of  H  =  H  .  The 

—  w  ’"e 

relation  between  0  and  0  is 


0  - 


(4) 


The  primes  in  Equations  (2)  and  (3)  indicate  properties  evalu¬ 
ated  at  the  local  reference  enthalpy,  given  by  the  expressions 
(Refs.  6  and  7) 


H'  *  0.23  He  +  0.19Hr+  0 . 58!  ^  ( laminar)  (5) 

H'  =*  0.36  He  +  0.19:Hr  +  0 . 45  Hw  ( turbulent)  (6) 


where 
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and 


(Pr  ■  )  1/2 

(laminar) 

(7) 

(Pr '  )  x^3 

(turbulent) 

(8) 

The  first  step  in  the  solution  of  Equation  (1)  consists  of 
determining  an  expression  for  the  right-hand  side,  that  is,  the 
Stanton  number.  To  do  this,  we  find  the  relation  existing  between 
0  and  CH  in  the  flat-plate  case,  and  then  assume  that  this 
relation  also  holds  for  the  body  of  revolution  in  a  pressure 
gradient . 

For  the  flat-plate  case  the  following  relations  are  known  to 

hold : 


C 


H 


where  for  laminar  flow 


(9) 


a  =  0.332 
m  =  0.5 

and  for  turbulent  flow  in  the  Reynolds  number  range  of  interest  for 
rocket  nozzle  application 


a  =  0.0296 


m  =  0.2 


Combining  Equations  (2)  and  (9) ,  there  results  for  the  flat-plate 
case 


p  u  (H  -  H  ) 
re  e  r  w 


a(p')1~n,(u')n; 

puV(Pr‘)2/3 
r  e  e 


(1C) 


Inverting  the  equation,  there  results 


m 


x  = 


m 

i— in 


(17) 


Combining  Equation  (17)  with  (10)  yields  an  expression  for  the  heat 
flux  term  that  does  not  contain  x  explicitly. 


p  u  (H  -  H  ) 
re  e  r  w 


=  C 


.  £l 


H 


(18) 


(Pr') 


It  is  now  assumed  that  Equation  (18) ,  which  was  derived  for 
the  case  of  the  flat  plate,  also  holds  for  the  body  of  revolution 
in  an  axial  pressure  gradient.  We  are  saying  that  the  energy 
thickness  and  the  heat  flux  are  related  in  the  same  manner  in  the 
two  cases.  Thus,  combining  Equation  (18)  with  (1),  there  results 
for  the  body  of  revolution: 
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Or, 


rewritten , 
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4“  rp  u  (H  -  H  ) 
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This  is  a  first-order  differential  equation  of  the  Bernoulli  type, 
and  can  be  solved  by  transforming  the  dependent  variable,  </>. 

Let 
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Then 


(1  -  m)p-m  +  -—A. 
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where  g  represents  the  bracketed  term  on  the  right-hand  side  of 
Equation  (20) . 

Equation  (22)  then  simplifies  to 

dP  +  (l-rrOrp^fH  -  H)  dZ  [rPeUe (Hr  "  Hw}]  P  dx  “  T^m  dx 
e  e  r  w 


(23) 


The  integrating  factor  for  this  equation  is 
/  (1-nOrp  u  \h  -  H  )  d  [rPeUe(Hr  "  Hw}] 


ln  -  Kv,  -  V]  " 


(24) 
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Hence ,  Equation  (23)  becomes 
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where  X  is  a  dummy  length  variable  along  the  surface.  Recalling 
that 
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_  i-m 
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Equation  (25)  becomes 
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And  now  Equation  (27) ,  in  conjunction  with  Equation  (18) ,  is 

used  to  calculate  the  wall  convective  heat  flux  or  Stanton  number 

on  or  in  an  axisymmetric  body  in  compressible  flow.  The  relations 

are  also  applicable  to  two-dimensional  flow  if  r  is  replaced  by 

unity  throughout.  The  integral  in  Equation  (27)  must  be  evaluated 

numerically,  but  it  is  a  simple  quadrature  and  can  be  handled 

quite  well  with  Simpson's  rule.  The  wall  temperature  must  be 

known  before  the  calculation  can  proceed,  since  the  reference 

properties  and  Hw  are  dependent  on  it.  Unfortunately,  the  wall 

temperature  is  often  precisely  what  we  are  trying  to  calculate,  so 

an  iterative  procedure  is  called  for.  A  trial  wall  temperature 

distribution  is  assumed,  CH'  is  calculated,  then  q^,  then  Tw 

(with  a  wall  conduction  program) ,  etc.  Normally  C  '  is  not 

n 

strongly  dependent  on  T^,  so  convergence  is  rapid. 

Another  quantity  required  to  initiate  the  calculation  is  0.^. 
If  no  other  information  is  available  for  estimating  0.  in  any 
given  application,  it  is  recommended  that  it  be  calculated  by 
equating  a  form  of  the  simple  Bartz  equation  and  Equation  (18)  at 
the  nozzle  entrance.  The  simple  Bartz  equation  form  is 


Equating  Equation  (28)  with  (18) ,  and  rearranging,  yields 
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2. 2. 1.2  Comparison  of  the  Ambrok-Rubesin-Mayer  results  with 
experimental  data 

The  data  used  to  check  the  ARM  method were  obtained  at  the 
Jet  Propulsion  Laboratory  and  are  reported  in  Reference  8.  The 
experimental  apparatus  used  is  shown  in  Figure  3 ,  taken  from  Ref¬ 
erence  8.  Air  was  heated  by  lean  combustion  of  methanol  in  a 
region  far  removed  from  the  nozzle  test  section  and  separated  by 
baffles,  screens,  and  a  long  duct.  The  nozzle  used  has  a  throat 
diameter  of  1.803  inches,  a  length  of  5.925  inches,  a  contraction 
area  ratio  7.75  to  1,  an  expansion  area  ratio  of  2.68  to  1,  a 
convergent  half-angle  of  30°,  and  a  divergent  half-angle  of  15°. 

It  was  fabricated  of  502-type  stainless  steel,  and  has  embedded 
thermocouples  for  measuring  local  heat  flux,  and  wall  static  pres¬ 
sure  taps.  The  heat  flux  is  obtained  from  the  temperature  drop  in 
the  stainless  wall  during  steady  state  water  jacket  cooling. 

The  particular  test  for  which  sufficient  data  was  reported 
to  permit  comparison  was  test  number  218.  In  this  test  the  stag¬ 
nation  pressure  was  75.2  psia,  and  the  stagnation  temperature  was 
1518°  R.  The  experimental  uncertainty  in  the  measured  heat  fluxes 
was  estimated  to  be  +11  percent  at  the  throat,  and  possibly  larger 
at  other  regions. 

The  nozzle  wall  axial  temperature  distribution  was  reduced 
from  a  temperature  profile  plot  given  in  Reference  8,  and  is  pre¬ 
sented  here  in  Figure  4.  The  measured  wall  pressure  distribution 
is  shown  in  Figure  5-  It  differs  significantly  from  that  calcu¬ 
lated  from  one-dimensional  flow  considerations,  indicating  signifi¬ 
cant  two-dimensional  effects  in  this  small  nozzle.  The  measured 
pressures,  and  the  corresponding  velocities  (assuming  isentropic 
expansion),  were  used  in  the  calculations. 

The  air  specific  heat  at  constant  pressure  varied  over 
the  temperature  range  of  interest,  so  a  curve-fit  expression  was 
used  to  relate  temperature  and  enthalpy. 
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H  -  14.6  +  0.2122T  +  16.93xl0-s  T2  -  2668.7  T 

where  H  is  given  in  Btu/lb  and  T  in  °R. 

Since  the  combustion  mixture  was  lean,  the  Sutherland  relation  for 
air  was  used  to  calculate  the  viscosity. 

M.  =  7.30x10  —  +  198  6 

where  |x  is  given  in  lb/sec-ft  and  T  in  °R. 

Since  the  boundary- layer  displacement  thickness  at  the 
* 

nozzle  entrance,  5^,  was  measured  during  the  JPL  tests,  this  value, 

rather  than  Equation  (29),  was  used  to  estimate  0..  The  measured 
•k  i 

value  of  6.  was  0.041  inch,  and  the  tables  of  Reference  9  lead 
1  *  _ 

to  an  estimated  value  of  (5  / ^  of  0.529.  Hence 

0.  =  0.00645  ft. 

The  corresponding  value  resulting  from  Equation  (29)  is  0.00182  ft. 
The  significance  of  this  discrepancy  will  be  discussed  after  the 
presentation  of  computed  results. 

Equation  (27)  was  solved  numerically  for  several  values 
of  0^,  including  the  "best  estimate"  of  0.00645  ft,  and  the 
results  are  shown  in  Figure  6.  The  effect  of  the  kink  in  the  pres¬ 
sure  distribution  at  x/L  =  0.7  can  clearly  be  seen.  Equation  (18) 
was  then  used  to  calculate  q^,  and  these  results  are  presented 
in  Figure  7,  also  for  several  values  of  0^.  Again,  the  effect  of 
the  pressure  kink  can  be  seen.  Especially  significant  is  the  strong 
effect  of  the  initial  energy  thickness,  0^. 

Figure  8  presents  the  comparison  of  t.he  calculated  q^ 

( 0^  =  0.00645  ft)  with  the  JPL  experimental  results.  Included  on 
the  figure  is  the  heat-flux  distribution  obtained  from  the  Bartz 
Equation,  (28) .  It  can  be  seen  that  the  Bartz  equation  results 


l 
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are  some  40  percent  above  the  ARM  results  and  the  experimental 
data.  Figures  9  and  10  show  the  calculated  enthalpy  convective 
heat-transfer  coefficient  and  Stanton  number,  respectively. 

It  can  be  concluded  from  Figure  8  that  for  the  experi¬ 
mental  conditions  considered  the  ARM  method  does  a  much  better 
job  of  predicting  the  data  than  does  the  Bartz  equation.  However, 
in  view  of  the  great  dependence  of  the  ARM  results  on  0k ,  as 
shown  by  Figure  7 ,  the  superiority  of  the  method  can  logically  be 
questioned  in  cases  where  0k  is  unknown.  As  discussed  earlier, 
when  Equation  (29)  is  used  to  estimate  0  for  the  case  at  hand, 

the  resulting  value  is  0.00182  feet  rather  than  0.00645  feet.  This 
leads  to  a  calculated  q^  about  25  percent  above  the  data  (see 
Figs .  7  and  8) . 

It  is  probable  that  Equation  (29)  will  do  a  better  job 
of  estimating  0^  in  cases  where  the  wall  is  not  as  highly  cooled 
(solid  rocket  nozzles)  ,  since  the  wall  cooling  increases  the  energy 
thickness;  however,  if  this  is  the  case,  the  Bartz  and  ARM  q^ 
results  would  also  be  closer  together.  Thus,  an  unqualified  recom¬ 
mendation  of  the  ARM  method  will  have  to  await  comparison  with 
experimental  data  taken  under  various  conditions  of  Tq/Tw-  This 
will  be  done  during  the  next  quarter. 

2.3  Surface  Chemical  Reactions  at  a  Graphite  Wall 

2.3.1  Method  of  analysis 

With  a  knowledge  of  the  transfer  coefficients,  we  seek  to 
generate  the  boundary  conditions  which  are  needed  for  the  solution 
of  the  temperature  history  of  the  wall  material.  Two  specific 
items  are  sought : 

(1)  The  net  heat  flux  to  the  surface. 

(2)  The  wall  erosion  rate. 

Both  these  quantities  are  in  general  found  as  a  function  of  wall 
surface  temperature,  and  thus  are  not  tied  down  for  any  specific 
application  until  the  resulting  relations  are  coupled  with  the 
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wall  temperature  history  calculation  discussed  in  Section  2.4.  To 
calculate  the  first  quantity  we  will  need  to  know  the  second  and 
the  enthalpy  of  the  gas  mixture  at  the  surface.  The  wall  erosion 
rate  enters  into  the  energy  balance  and  the  gaseous  products  modify 
the  heat-transfer  coefficient.  To  calculate  these  we  turn  to  an 
analysis  of  the  chemical  reactions  occurring  at  the  exposed  surface. 

The  exposed  surface  of  a  nozzle  liner  material  is  the  scene 
of  intense  chemical  activity.  For  example,  for  a  graphite  wall  in 
a  typical  H,  C,  0,  N,  Cl  exhaust  there  are  about  a  dozen  molecular 
species  which  appear  in  significant  quantities  during  the  course  of 
a  firing.  Hence,  there  is  a  large  number  of  chemical  reactions 
which  can  (and  undoubtedly  do)  occur.  The  task  of  analyzing  these 
reactions  is  further  complicated  by  the  fact  that  not  all  of  them 
achieve  equilibrium  under  the  conditions  of  interest. 

To  illustrate  the  implications  of  "non-equilibrium"  in  a 
flow  system  which , nevertheless ,  can  be  essentially  steady  state, 
consider  the  control  volume  shown  in  the  following  sketch: 


Gaseous 

boundary 

layer 


Surface 

material 
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The  lower  boundary  of  the  control  volume  is  taken  right  at  the 
surface  of  the  material  exposed  to  the  exhaust  stream.  The  con¬ 
trol  volume  is  presumed  to  be  only  infinitesimally  thick.  Assum¬ 
ing  steady  state,  that  is,  no  mass  storage  in  the  control  volume, 
a  mass  balance  for  any  molecular  species  can  be  written  as  follows 


(pv)  K.  +  p  u  C.,  (K.  -  K.  )  =  m  K.  +  w, 

r  w  iw  re  e  M.  iw  le  v  iv  i 


(30) 


where 


(pv), 


net  convection  mass  velocity  normal  to  the  surface 


K. 


xw 


mass  fraction  of  species  i  in  the  surface  gas 
mixture 


'M. 

1 


mass-transfer  coefficient  of  species  i  through 
the  mixture 


K 


ie 


m 


mass  fraction  of  species  i  on  the  free  stream 
erosion  rate  of  virgin  wall  material 


K. 

xv 


mass  fraction  of  species  i  in  the  wall  material 


w . 

l 


net  rate  of  production  of  species  i  due  to 
chemical  reactions 


The  use  of  the  convective  mass-transfer  term  in  Equation  (30) ,  and 
in  particular  the  use  of  the  difference  in  molecular  (rather  than 
atomic)  concentration  as  the  potential  for  mass  transfer  merits 
some  discussion.  Strictly  speaking,  this  expression  should  apply 
only  to  a  nonreacting  boundary  layer,  and  it  is  not  unreasonable 
here  to  impose  the  assumption  that  the  boundary  layer  is  frozen. 

It  may  not,  however,  be  necessary  to  make  this  assumption. 

Kulgein  (Ref.  27)  measured  mass-transfer  rates  for  methane  with 
burning  in  the  boundary  layer.  He  found  that  the  wall  to  stream 
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partial-pressure  difference  correlated  his  results  both  with  and 
without  reaction.  Thus,  the  possibility  exists  that  a  less 
restrictive  assumption  than  that  of  a  frozen  boundary  layer  will 
still  permit  Equation  (30)  to  be  valid,  and  this  will  be  explored 
further  during  the  coming  quarter. 

The  chemical  production  term  in  Equation  (30)  can  be  elimi¬ 
nated  by  converting  to  an  elemental  mass  balance;  this  is  done  by 
summing  the  mass- balance  equations  (appropriately  weighted)  of  all 
the  molecular  species  which  contain  a  given  element.  The  result  is 
a  set  of  equations  relating  the  elemental  mass  fractions  at  the  wall, 
K^w,  to  the  composition  of  the  free  stream  with  the  erosion  rate 
and  mass-transfer  coefficients  as  parameters.1  However,  in  order 
to  relate  erosion  rate  to  temperature  and  pressure,  and  to  enable 
evaluation  of  the  wall-gas  enthalpy,  we  need  to  know  the  molecular 
composition,  that  is,  the  K.  .  To  do  this  we  must  write  Equa¬ 
tion  (30)  for  all  the  species  of  interest  and  solve  the  set  of 
equations  simultaneously. 

The  nature  of  the  chemical  production  term  is  best  illus¬ 
trated  by  a  specific  example.  For  the  reaction 


2C  +  H, 


C2H2 


rate  data  indicate  that  the  production  of  C2H2  at  surface 

can  be  written  as 


C2H2 


kfPH. 


production 


where 


reaction  rate  constant 


PH  partial  pressure  of  H2  at  the  surface 


This  is  discussed  in  detail  in  Reference  1. 
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If  we  assume  that  the  reverse  reaction  can  be  written  as 


consumption 


then  the  definition  of  equilibrium  provides  us  with 


(33) 


where  is  the  equilibrium  constant  for  the  reaction  given  by 

Equation  (31) ,  and  is,  of  course,  a  function  of  temperature.  Hence, 
the  net  production  rate  is  given  by 


For  a  closed  system,  that  is,  a  system  consisting  of  a 
mixture  in  a  box,  regardless  of  the  value  of  kf,  equilibrium  is 
obtained,  if  one  is  patient  enough,  and  w  vanishes.  The  control 
volume  shown  in  the  previous  sketch  is  an  open  system  in  the  sense 
that  mass  crosses  the  boundaries.  For  such  an  open  system  equi¬ 
librium  is  never  attained,  it  is  only  approached.  Further,  the 
approach  to  equilibrium  occurs  when  the  reaction  rates  are  large, 
that  is,  when  kf  — *■  <».  As  a  result  of  this  w  remains  finite 
as  equilibrium  is  approached  in  the  open  system.  In  fact,  w 
vanishes  only  for  k^  — +  0 ,  that  is,  when  no  reactions  occur. 
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To  illustrate  the  behavior  of  Equation  (30)  as  —• ►  <», 

substitute  Equation  (34)  into  Equation  (30)  and  rearrange  somewhat. 
We  have  then2 


'  C_H 


2  2 


K 


PeUeCM 


(35) 

As  becomes  very  large  the  right-hand  side  of  Equation  (35) 

approaches  zero,  which  results  in 


■  C_H 


2  2 


K 


(36) 


This  is  exactly  the  equilibrium  relation  between  the  partial  pres¬ 
sures  of  the  species  involved  in  the  reaction  of  Equation  (31) . 

Then  as  the  reaction  rates  become  large,  the  mass-balance  equations, 
Equation  (3C) ,  reduce  to  the  equilibrium  relations  among  the  product 
species.  Notice  that  this  same  result  will  obtain  even  if  the 
system  behaves  non-ideally,  that  is,  the  order  of  the  reaction  is 
not  givWi  by  the  stoichiometric  coefficient. 

If  sufficient  rate  data  are  available  the  system  of  equa¬ 
tions  can  be  set  up  and  solved.  Let  us  review  the  steps  involved 
in  such  a  calculation. 

(a)  Specify  problem  parameters: 


w' 


■  w1 


p  u  C„  , 
re  e  M.  ■ 

i 


K.  , 
le 


K .  , 

le  * 


kf  i 


2Note  that 


0. 
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(b)  Assume  (pv)^  and  calculate  Kiw  using  the  elemental 
mass  balance  equations  resulting  from  the  weighted  sums  of  Equa¬ 
tion  (30) - 3 

(c)  Calculate  the  K^w  using  the  set  of  equations  resulting 
from  Equation  (30) ,  and  the  expression  for  conservation  of  elemental 
species , 


K. 


1W 


I 


a 

ij 


(37) 


where  is  the  mass  fraction  of  element  i  in  molecule  j . 

Steps  (b)  and  (c)  are  repeated  in  an  iterative  manner  until 
the  system  closes,  that  is,  until  all  of  the  equations  resulting 
from  Equation  (30)  and  Equation  (37),  are  satisfied.  The  result  is 
a  unique  value  for  (pv)^  and  since  the  are  known  the  wall- 

gas  enthalpy  can  be  calculated. 

We  have  then  in  principal  solved  the  system  and  generated 
the  desired  boundary  conditions  for  the  wall  temperature  calcu¬ 
lation.  There  are,  however,  a  number  of  difficulties  which  merit 
discussion.  They  are: 

(1)  Rate  data  for  many  of  the  reactions  of  interest  are 
either  non-existent  or  are  singularly  unreliable. 

(2)  Even  if  perfectly  reliable  data  were  available  the 
problem  involves  many  parameters  and  a  concise  general  description 
of  the  solutions  is  not  possible. 

(3)  The  mass-transfer  coefficients  in  a  reacting  system 
are  not  well-established  experimentally  although  recent  work  has 
shed  light  on  this  subject. 

Items  (1)  and  (2)  will  be  discussed  here,  item  (3)  will  be  con¬ 
sidered  in  a  later  report. 


3Details  of  this  step  are  presented  in  Reference  1. 
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2.3.2  Reaction  kinetics 

The  determination  of  reaction  rates  for  heterogeneous 
reactions,  particularly  at  high  temperatures,  is  at  best  a  diffi¬ 
cult  task.  As  temperature  is  increased,  the  mixture  approaches 
an  equilibrium  composition  at  the  reacting  surface  and  becomes 
rate  limited  by  the  transport  of  reactant  to  the  surface  (diffusion 
control) .  As  can  be  seen  from  Equation  (35)  the  degree  of  approach 
to  equilibrium  depends  in  part  on  the  magnitude  of  the  mass-transfer 
coefficient.  Since  the  mass-transfer  coefficients  obtained  in 
rocket  nozzles  are  very  high,  the  onset  of  diffusion  control  is 
delayed  and  the  desired  data  are  in  a  range  unobtainable  by  most 
conventional  rate  experiments.  Further  difficulties  in  determi¬ 
nation  of  reaction  rates  are  posed  by  the  fact  that  the  physical 
character  of  the  reacting  surface  can  have  a  considerable  influence 
on  the  reaction  rates  and  that  the  accurate  measurement  of  surface 
temperature  is  difficult , particularly  at  high  temperatures. 

The  combination  of  these  factors  causes  what  data  that  do 
exist  to  be  suspect  and  often  contradictory.  An  illuminating 
discussion  of  existing  rate  data  for  the  reaction 

2C*  +  02  ^  2C0  (38) 

is  given  by  Scala,  Reference  10.  In  summarizing  the  data  of 
numerous  workers,  Scala  points  out  that  the  order  of  the  reaction 
is  variously  reported  as  lying  between  zero  and  one,  the  activa¬ 
tion  energy  resulting  from  assuming  a  standard  Arrhenius  form 
varies  from  8  to  60  Kcal/mole  and  the  pre-exponential  factor  varies  — 
over  several  orders  of  magnitude. 

An  insidious  complication  in  heterogeneous  reactions  is  the 
fact  that  in  order  to  react,  a  molecule  must  get  to  the  surface. 

In  a  mixture  of  gases,  such  as  a  rocket  exhaust,  there  is  compe¬ 
tition  among  the  various  reacting  species  and  indeed  among  molecules 
of  the  same  species  for  available  surface  sites.  In  fact,  even  an 
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inert  diluent  can,  if  strongly  adsorbed,  inhibit  surface  reactions. 
The  net  result  is  that  reaction  rates  in  a  mixture  of  gases  can  be 
slower  than  rates  measured  singly  or  at  very  low  pressures. 

Some  methods  of  correcting  for  this  interaction  effect  have 
been  summarized  by  Trapnell  (Ref.  11),  both  for  mobile  and  immobile 
layers.  For  example,  assuming  a  mobile  layer,  that  is,  one  in  which 
the  molecules  may  move  easily  from  site  to  site,  the  effect  results 
in  a  change  in  activation  energy  of  the  reaction  given  by 


AE  =  asV 


1  - 


l  -  2n 


[l  -  4(1  -  Oti  (1  -  rj)]  1 


71 


(39) 


where 


s 

V 

T] 

C 

a 


number  of  neighboring  sites 

repulsive  potential 

fraction  of  sites  covered 

-V/RT 

e 

multiplying  factor;  unity  if  -the  reaction  requires  two 
surface  sites,  one-half  if  only  one  is  needed 


As  an  upper  limit  T]  — *■  1 ,  and  with  s  =  4  and  a  =  1,  (which  would 
represent  a  two  site  reaction  on  graphite) ,  Equation  (39)  reduces 
to 


AE  =  8V 


As  an  estimate  of  this  effect,  Roberts  (Ref.  12)  gives  AE  as 
27  Kcal  for  the  adsortpion  of  hydrogen  molecules  on  tungsten. 

A  more  severe  interference  with  reaction  occurs  if  the 
surface  is  wholly  covered  by  species  which  act  as  poisons  and  must 
be  removed  before  reaction  can  occur.  Following  Glasstone, 

Laidler,  and  Eyring  (Ref.  13) ,  one  can  express  the  reaction  rate  in 
such  a  situation  as 
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(w)  f  = 


ZU 


1=J 


(40) 


where 

(w) ^  a  uni-directional  or  forward  rate 


partial  pressures  of  poisons 
Pj  partial  pressure  of  reactant 


The  rate  constant  kf  involves  a  pre-exponential  factor  that  may 
be  estimated  theoretically,  but  the  activation  energy  requires 
experimental  determination. 

These  considerations  indicate  further  the  difficulty  in 
obtaining  heterogeneous  rate  data.  One  never  has  pure  reactant  at 
the  surface,  the  inhibiting  effect  of  the  product  species  must 
always  be  considered,  and  dissociation  at  high  temperatures  may 
add  other  molecular  species  to  the  mix  at  the  surface. 


2.3.3  Surface  recession  of  a  graphite  wall 

A  series  of  calculations  was  performed  for  a  pure  graphite 
wall  in  a  typical  exhaust  stream  excluding,  however,  all  Al^O.^. 
The  purpose  of  these  calculations  was  two-fold,  first  to  assess 
the  contribution  of  the  various  surface  reactions  to  the  erosion 
of  the  wall  by  chemical  reactions,  and  second,  to  evaluate  the 
importance  of  the  various  parameters  affecting  the  calculation. 

In  the  absence  of  alumina  particles  there  are  three  prin¬ 
cipal  combustion  products  which  attack  the  wall  chemically.  The 
basic  reactions  are  (determined  from  equilibrium  chemistry  calcu¬ 
lations)  : 

C*  +  co2  ^ 


2C0 


(41) 
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The  same  set  of  possibilities  were  considered  for  reaction 
(43) .  The  rate  data  were  extrapolated  from  measurements  made  by 
the  Union  Carbide  Research  Institute  (Refs.  16  and  17).  An  Arrhenius 
plot  of  these  data  is  shown  as  Figure  11.  The  equation  for  the 
forward  reaction  rate  constant  was  taken  as 


kf  =  4.525x10s  e 


lb/ ft 2-sec-Atmos 


(47) 


where  T  is  in  °K.  It  should  be  noted  that  the  relia¬ 
bility  of  the  surface  temperature  measurement  (which  was  a  "bright¬ 
ness"  temperature  measured  by  a  pyrometer)  is  difficult  to  assess. 

For  each  of  the  three  possibilities  for  reactions  (42)  and 
(43) ,  comparison  calculations  were  made  to  obtain  the  dimensionless 
wall  recession  rate 


B' 


(Pv)w 

P  u  c», 
re  e  M 


(48) 


as  a  function  of  wall  temperature  for  the  following  conditions 


p  =10  atmospheres 

PeugCM  =  1  (and  is  the  same  for  all  species) 


0.1178 


C2H2. 


The  rate  constants  were  calculated  from  Equations  (45) ,  (46) ,  and 

(47).  The  effect  of  the  reverse  reactions  was  ignored,  however, 
the  forward  reaction  was  constrained  so  as  not  to  proceed  beyond 
equilibrium.  Under  this  assumption  equilibrium  is  approached  more 
rapidly  than  it  should  be.  An  estimate  for  the  effect  of 
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interaction  of  competing  species  was  made  by  increasing  the  acti¬ 
vation  energies  by  24  Kcal/mole  (a  change  of  this  order  would  be 
indicated  for  an  immobile  layer) .  Because  there  is  some  hydro¬ 
gen  dissociation,  particularly  at  the  high  end  of  the  temperature 
range  considered,  the  possibility  exists  that  the  reaction 

2C *  +  2H  ^  C2H2  (49) 

proceeds  more  rapidly  than  reaction  (43) .  Accordingly,  the  rate 
constant  for  reaction  (49)  was  arbitrarily  taken  as  10  times  that 
given  by  Equation  (47) . 

In  addition  to  C2H25  the  f°ll°win9  hydrocarbon  molecules 
were  considered;  CH4 ,  CH2 ,  C2H,  C3H,  C4H,  C6H,  C4H2 ,  and  C^. 

The  methane  was  assumed  to  be  in  equilibrium  with  the  graphite  and 
the  hydrogen  present  at  the  surface  (very  little  methane  remains 
above  2500  K  anyway) ,  while  the  other  C-H  molecules  were  assumed 
to  be  in  equilibrium  with  graphite  and  C2H2 ,  thus  tying  their 
rate  of  production  to  that  of  acetylene.  The  latter  assumption 
is  somewhat  artificial  and  highly  questionable.  Incomplete  infor¬ 
mation  reported  in  Reference  17  indicates  that  the  molecular  frag¬ 
ments  considered  may  form  much  more  slowly  than  the  rate  required 
to  equilibrate  with  C2H2. 

The  remaining  molecular  species  were  allowed  to  distribute 
according  to  the  dictates  of  chemical  equilibrium.  The  propellant 
composition  and  a  complete  list  of  product  species  considered  are 
shown  in  Table  I.  These  calculations  were  made  by  a  digital  com¬ 
puter  program  prepared  for  Vidya ' s  IBM  1620. 

Figure  12  is  a  plot  of  B'  versus  temperature  for  some  of 
the  possible  alternative  situations.  It  is  readily  apparent  that 
rate  assumptions  regarding  just  these  two  reactions  can  profoundly 
affect  the  prediction  of  wall  chemical  erosion  rate.  At  3000°  k , 
for  instance,  the  assumption  that  neither  Equation  (42)  nor  (43) 
occur,  reduces  B'  to  1/5  of  the  value  predicted  by  equilibrium. 
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If  the  rate  data  given  here  are  realistic,  B1  is  only  0.4  to  0.5 
of  the  value  predicted  by  equilibrium.  Thus,  it  is  clear  that  there 
is  a  real  need  for  more  complete  and  carefully  acquired  fundamental 
rate  data. 

Assuming  kinetic  control  of  reactions  (42)  and  (43) ,  the 
rate  constants  given  by  Equations  (45) ,  (46) ,  and  (47) ,  together 

with  the  same  assumptions  made  above,  a  series  of  calculations  was 
made  to  study  the  effects  of  temperature,  pressure,  and  mass-transfer 
coefficient.  These  results  are  presented  as  Table  II.  The  evalu¬ 
ation  of  these  results  is  still  in  progress,  but  Figure  13  illus¬ 
trates  the  magnitude  of  the  effects  which  may  be  encountered.  The 
range  of  variables  considered  was  as  follows : 

5.0  <  p  <  50.0  (atmospheres) 

1500  <  T  <  3600  (°K) 

0 . 1  <  p  u  C  <  2.0  ( lbs/ft 2-sec) 

—  '  e  e  M  — 


At  any  given  temperature  both  the  system  pressure  and  the  mass- 
transfer  coefficient  can  significantly  affect  B'.  For  example, 
at  3000°  K  increasing  pressure  from  5  to  50  atmospheres  increases 
B'  by  about  1/3.  At  the  same  temperature  an  increase  in  PeueCM 
of  from  0.1  to  1  decreases  B’  by  a  factor  of  2.  It  should  be 
noted,  however,  that  since  B’  contains  the  factor  p  u  C„  in 
the  denominator  this  really  represents  an  increase  in  wall  erosion 
by  a  factor  of  5. 

An  effect  not  illustrated  by  these  calculations  is  the  sen¬ 


sitivity  of  B' 


to 


(>2°) 


If  the  rate  of  reaction  (42)  is 


zero,  increasing  the  mass  fraction  of  water  in  the  free  stream 
from  0.1178  to  0.13  (as  can  occur  with  a  shift  in  stream  tempera¬ 
ture)  decreases  B‘  from  0.027  to  0.016.  This  redaction  occurs 
because  both  the  oxygen  and  the  nydrogen  present  in  the  exhaust 


gas  as  water  vapor  are  prevented  from  reacting  with  the  wall.  This 
effect  is,  of  course,  less  pronounced  as  reaction  (42)  nears  com¬ 
pletion,  since  the  amount  of  remaining  water  vapor  in  equilibrium 
at  the  wall  is  small.  It  is  fortunate  that  if. the  exhaust  stream 


is  in  equilibrium  (or  near  it)  the  value  of 
zero,  because  a  similar  dependence  of  b' 


^”2^2 


is  nearly 


on  the  acetylene  fraction 


exists . 


Figure  13  together  with  the  last  remarks  indicate  the  diffi¬ 
culty  in  providing  the  wall  boundary  conditions  for  the  conduction 
solution,  in  that  pressure,  the  mass-transfer  coefficient,  and  the 


amount  of  free  stream  H 2 O  will  vary  both  with  position  in  the 
nozzle  and  with  time.  Two  approaches  are  being  explored: 

(1)  Include  the  calculation  of  the  surface  chemical  reac¬ 
tions  as  a  subroutine  to  the  transient  conduction  program. 

(2)  Develop  approximate  correlations  so  that  a  reasonable 
table  look-up  routine  can  be  used. 

2.3.4  Graphite  ablation  in  Mixture  4 

As  discussed  in  detail  in  Section  3,  the  test  firings  of  a 
graphite  nozzle  using  an  He,  N2 ,  02  mixture,  Mixture  4,  have 
already  begun.  If  we  accept  the  premise  of  equal  mass-transfer 
coefficients  for  all  species  and  assume  that  the  reaction  (38) 
achieves  equilibrium  under  the  test  conditions,  the  wall  erosion 
rate  assumes  a  particularly  simple  form. 

In  the  approximate  wall  temperature  range  2000°  K  to  3300°  K, 
equilibrium  predicts  that  the  only  significant  molecular  species 
at  the  wall  are  He,  N2 ,  and  CO  while  the  main  stream  contains  only 
He,  N2,  02 ,  and  O.  Then  the  oxygen  balance  can  be  written  as6 


and  the  carbon  balance  as 

sSee  Reference  1  for  a  derivation  of  this  form 


(50) 
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CO 


B' 


C  J  1  +  B' 

w 


But  since 


and 


GO  w  ■  a  CK0» 


.  16  Gk  "0 

28  \^KCO/ 


w 


4  V  O 


w 


B* 


w 


1  +  B' 


Substituting  Equation  (52)  into  (50)  and  solving  yields 


B*  -  4  KC 


(51) 


(52) 


(53) 


so  that  in  order  to  estimate  the  wall  recession  rate  we  need  only  to 
know  the  composition  of  the  gas  mixture  and  the  local  mass-transfer 
coefficient.  It  should  be  emphasized  that  this  expression  is  for 
chemical  equilibrium  with  a  graphite  wall  and  applies  only  for  Mix¬ 
ture  4  at  wall  temperatures  between  2000°  K  and  3300°  K.  Its  sim¬ 
plicity  is  due  to  the  fact  that  only  the  oxygen  is  chemically  reactive 
in  this  temperature  range.  The  erosion  rate  for  chemical  equilibrium 
is  then  essentially  independent  of  wall  temperature  for  this  special 
case . 


2.3.5  Future  effort 

An  immediate  task  is  to  make  modifications  to  the  program 
for  calculating  the  chemical  erosion  rate  of  graphite.  The  effect 
of  reverse  reaction  will  be  included.  It  is  also  intended  to  relax 
the  equilibrium  assumption  for  reaction  (41) .  Further  consideration 
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will  be  given  to  the  phenomena  of  interaction  and  poisoning  on  the 
surface.  In  addition,  the  problem  of  providing  the  results  of  the 
chemical-erosion  calculation  in  an  optimum  form  for  use  in  the  con¬ 
duction  program  will  be  pursued. 

Considerable  effort  will  be  placed  on  an  analysis  of  the 
results  of  the  chemical  erosion  tests  being  conducted  in  the  Vidya 
plasma  facility.  It  may  be  possible  to  obtain  reaction  rate  data 
directly  from  these  tests  and  they  v<ill  certainly  provide  consider¬ 
able  guidance  for  the  selection  of  data. 

The  next  major  task  of  the  program,  the  analysis  of  char 
producing  ablators,  will  begin  during  the  next  quarter. 

2.4  Internal  Material  Thermal  Behavior 

As  indicated  in  Figure  2,  the  determination  of  the  internal 
material  thermal  behavior  is  an  integral  part  of  the  overall  wall- 
recession  calculation,  in  that  it  provides  the  wall  surface  tempera¬ 
ture  input  for  the  surface  chemical  reaction  evaluation.  This  wall 
response  is  calculated  with  a  transient  conduction  computer  program 
which  was  discussed  in  some  detail  in  the  First  Quarterly  Progress 
Report.  Further  work  on  this  computer  program  is  presented  in  this 
section. 

2.4.1  The  heat  equation  with  a  moving  boundary 

To  facilitate  the  numerical  solution  of  the  heat  equation  when 
one  boundary  surface  is  receding  due  to  ablation,  the  equation  for 
conduction  of  heat,  in  cylindrical  coordinates,  can  be  transformed 
as  follows.  If  we  set 

x  =  r  -  rg (z , 0) 
z 1  =  z 
0'  -  6 


the  heat  equation  becomes 


‘(54) 


The  motives  for  the  selection  of  this  form  are  discussed  in  detail 
in  Reference  1. 

For  computation,  the  various  derivatives  appearing  in 
Equation  (54)  are  expressed  in  finite  difference  form.  The  nodal 
network  is  illustrated  in  the  following  sketch. 
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Since  the  problem  is  axisymmetric ,  the  wall  section  shown 
is  actually  a  section  of  a  body  of  revolution.  The  coordinate 
system  used  poses  some  restrictions  on  the  placement  of  nodes.  A 
row  of  nodes  always  appears  on  the  surface.  Az (m)  can  be  unequal 
to  Az(m-l)  but  Az(m)  is  a  constant  for  all  n.  Similarly  Ax(n) 
need  not  equal  Ax(n-l)  but  Ax(n)  is  a  constant  for  all  n. 

The  particular  finite  difference  forms  of  the  various 
derivatives  in  Equation  (54)  are  lengthy  and  are  not  included  in 
this  report. 

As  discussed  in  Reference  1,  the  formulation  of  the  differ¬ 
ence  equation  representing  Equation  (54)  can  be  either  explicit  or 
implicit.  The  program  for  the  solution  of  Equation  (54)  utilizing 
the  explicit  formulation  has  been  completed  and  checked  out. 

2.4.2  Ablating  surface  boundary  condition 

The  temperatures  for  nodes  not  located  on  (or  near)  a  bound¬ 
ary  are  calculated  from  the  difference  equation  representing  Equa¬ 
tion  (54) .  Special  considerations  are  required  for  those  points 
located  at  the  ablation  surface  and  there  is  some  latitude  in  the 
method  used  to  calculate  surface  temperatures.  Two  different 
methods  can  be  used:  (a)  heat-flux  balance  at  the  surface, 

(b)  energy  quantity  balance  on  a  finite  region  surrounding  a  sur¬ 
face  node. 

To  develop  some  quantitative  experience  for  the  relative 
merits  of  the  two  methods ,  the  one-dimensional  heat  conduction 
equation  was  programed  in  finite  difference  form  utilizing  a  simple 
convective  (constant  h)  boundary  condition  and  constant  material 
properties.  A  further  objective  in  constructing  the  simple  program 
was  to  permit  an  evaluation  of  the  size  of  the  finite  difference 
net  required  to  provide  a  good  approximation  to  the  differential 
equation . 
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Referring  to  the  representation  of  the  one  dimensional 
geometry  shown  below, 
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'convection 

the  boundary  condition  given  by  the  heat-flux  balance  takes  the 
form 


qconvection  h^Ta  Ti^  k 

If  the  temperature  gradient  at  the  surface  is  obtained  from  Taylor's 
series  expansions  from  which  the  second  derivative  terms  have  been 
eliminated,  one  obtains  the  finite  difference  form: 


h  (T  -  T  ) 

u  X 


(56) 


which  may  be  solved  for  T  directly. 

The  alternative  form  results  from  an  energy  balance  on  the 
dashed  box  shown  in  the  above  sketch.  We  have  then 

Ax  St  f  dT  \ 

2  c)0  ^convection  +  k\^x  ^  (57) 

xc°  2 

Again,  if  one  obtains  (dT/dx)  from  Taylor's  series  expansions,  and 
eliminates  the  second  order  terms,  there  results 


A0 


T 

_ l 


h  (T. 


-  V 


(T. 


+  k 


"  V 


Ax 


(58) 
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where  the  prime  denotes  temperature  at  the  end  of  the  time  interval, 
AS. 

Since  the  form  given  by  Equation  (56)  does  not  contain  a 
term  relating  to  energy  storage  it  is  a  relationship  that  applies 
at  any  instant  in  time.  Then,  after  the  calculation  of  the  "new" 
temperatures  for  all  the  interior  points.  Equation  (56)  can  be 
written  in  terms  of  the  end  of  time  interval  values,  T  ' ,  T  1 , 
and  T  1  and  solved  for  T  1  directly. 

3  1 

For  the  more  general  convective- boundary  condition,  the  net 
convection  heat-flux  is  written  in  terms  of  enthalpy  driving  poten¬ 
tials  and  contains  a  correction  due  to  surface  off-gassing.  In 
this  case,  an  iterative  procedure  will  be  required  to  solve  for  T  1 
and  it  becomes  sensible  to  consider  the  use  of  T^  on  the  left-hand 
side  of  Equation  (55)  while  retaining  T on  the  right.  This 
simplification,  however,  again  raises  the  problem  of  "stability" 
of  the  finite  difference  form.  This  is  best  illustrated  by  carrying 
out  this  simplification  as  follows.  The  expression  for  T  1 
becomes : 


(59) 


If  the  explicit  expression  for  T  '  is  substituted  into  Equa¬ 
tion  (59) ,  we  have 


(59  a) 


The  general  condition  of  stability  requires  that  the  coefficient  on 
T  appearing  in  Equation  (59a)  must  be  positive  (see  Appendix  A) . 
Thus,  the  inequality 
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aAg  k  1 

Ax3  h  >  2 


(60) 


must  be  satisfied  to  yield  a  meaningful  solution  for  surface  tempera¬ 
ture.  It  is  further  evident  that  if  the  convection  heat-flux  term 
is  written  in  terms  of  T  '  ,  Equation  (56)  is  unconditionally  stable 
for  all  values  of  A 6  and  Ax  because  the  term  containing  A 6 
and  T  which  enters  via  T  ’  is  positive.  Trial  solutions  in 
which  the  condition  given  by  (60)  was  violated  were  attempted  and 
did,  in  fact,  prove  troublesome. 

A  similar  decision  with  regard  to  the  application  of  Equa¬ 
tion  (58)  must  be  made.  If  the  formulation  is  explicit  (that  is, 

T i 1  appears  only  on  the  left-hand  side  of  (58)) ,  then  we  have  the 
familiar  condition  that 


QA  6 
Ax2 


h _  1 

(k/Ak)  <  2 


(61) 


while  the  implicit  expression  is,  again,  unconditionally  stable. 

Figure  14  is  a  comparison  of  the  calculated  surface  tempera¬ 
tures  resulting  from  the  finite  difference  scheme  for  each  of  the 
two  boundary  conditions  (56)  and  (58),  (using,  in  both  cases,  the 
stable  forms)  with  the  values  given  by  the  analytical  solution. 

Both  methods  give  results  somewhat  in  error  near  the  beginning  of 
the  problem  and  both  seem  to  approach  the  correct  solution  at  about 
the  same  time  (for  comparable  difference  mesh  size).  An  examination 
of  Equation  (56)  reveals  the  fact  that  at  the  end  of  the  first  time- 
step  the  same  result  is  obtained  regardless  of  the  value  of  A 6. 
Hence,  reducing  the  size  of  A Q  may  make  the  solution  worse  for 
short  times.  Equation  (58)  does  not  have  this  disadvantage  in  that 
reducing  A 6  always  brings  one  closer  to  the  correct  solution. 

On  this  basis,  Equation  (58)  has  been  used  in  setting  up  the  axi- 
symmetric  program,  but  the  situation  is  not  completely  settled. 
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In  order  to  derive  the  form  equivalent  to  Equation  (58)  for 
the  transformed  axisymmetric  coordinate  system,  we  must  include 
the  extra  terms  resulting  in  the  skewed  and  moving  system.  An 
energy  balance  on  the  surface  element  shown  in  the  sketch  on  page  32, 
yields  (assuming  Ax  and  Az  to  be  differentials) : 


where  the  terms  inside  the  braces  are  understood  to  be  evaluated  at 
x  =  dx/2  rather  than  at  the  surface. 

2.4.3  Program  verification 

Equation  (54)  together  with  the  boundary  condition  expressed 
by  Equation  (62)  have  been  orogramed  in  an  explicit  difference 
formulation.  To  verify  the  internal  logic,  this  program  was  used 
to  solve  the  problem  of  the  hollow  cylinder  (in  the  absence  of 
ablation)  having  constant  properties  and  a  simple  convective  bound¬ 
ary  condition.  The  results  of  the  numerical  solution  are  compared 
with  the  analytical  solution,  taken  from  Schneider  (Ref.  18),  in 
Figure  15.  As  can  be  seen,  the  comparison  is  very  favorable.  The 
problem  parameters  were: 

Ax  =  0.00833  feet 

A  =  0.1  sec 

r  =  0.125  feet 
s 


0.20833  feet 
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h  =  0.36  Btu/ft2-sec-^R 

k  =  0.0075  Btu/ft-sec-°P.  'l  typical  of 

pc  =  59.8  Btu/f t3- °R  J  graphite 

After  one  second  of  problem  time,  the  surface  temperatures  are 
within  1  percent  of  the  correct  values. 

The  next  step  in  program  checkout  is  the  construction  of 
an  ablating  wall  solution  from  the  non-ablating  wall  solution. 

This  is  done  as  follows.  Once  the  complete  solution  of  the  non¬ 
moving  wall  problem  is  at  hand,  any  arbitrary  surface  recession 
rate  can  be  specified.  Then,  since  the  position  of  the  surface  as 
a  function  of  time  is  known,  the  temperature  of  the  surface  can  be 
specified  as  that  temperature  which  existed  at  that  position  at 
that  time  in  the  non-ablating  body.  As  a  result  of  this,  the  portion 
of  the  body  which  remains  must  respond  exactly  as  if  no  ablation 
had  occurred  and,  for  example,  the  back  wall  temperatures  must  be 
identical  in  the  two  cases.  This  comparison  test  is  still  in 
progress . 

2.4.4  Future  effort 

The  most  immediate  task  is  the  ablating  wall  comparison  just 
discussed.  The  addition  of  the  more  general  boundary  conditions 
is  the  next  following  task.  These  two  items  would  complete  the 
transient  temperature  package  in  its  explicit  formulation. 

Because  of  the  mixed  derivative  terms  appearing  in  the 
transformed  heat  equation,  Equation  (54),  the  finite  difference 
expression  for  the  temperature  at  any  nodal  point  m,  n  involves 
the  array  of  temperatures  existing  at  the  point  m,n  and  the 
eight  surrounding  points.  In  preparing  the  coefficient  matrix  for 
the  implicit  solution,  the  existence  of  nine  non-zero  coefficients 
in  each  row  results  in  some  additional  complexity,  for  example,  the 
dominant  terms  must  be  located  and  moved  to  the  diagonal.  Efforts 
to  overcome  these  difficulties  look  hopeful  and  the  implicit  formu¬ 
lation  will  be  continued. 
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The  internal  behavior  of  a  charring  ablator  is  character¬ 
ized  by  two  moving  boundaries.  The  analysis  of  these  types  of 
materials  will  be  started  in  the  next  quarter. 

3.  EXPERIMENTAL  STUDIES 

3.1  Introduction 

The  choice  of  a  rocket-nozzle  liner  material  requires  a 
knowledge  of  the  behavior  of  appropriate  materials  in  the  combus¬ 
tion  product  environment-  This  knowledge  must  be  supplied  by 
both  an  analytic  treatment  of  the  materials  response  and  experi¬ 
mental  materials  testing.  As  indicated  in  Figure  2,  the  analytic 
treatment  requires  appropriate  input  of  basic  data  for  a  success¬ 
ful  and  accurate  solution,  and  some  of  this  input  must  be  deter¬ 
mined  experimentally.  Hence,  in  the  experimental  phase  of  this 
program,  several  areas  in  which  data  are  currently  lacking  are 
being  investigated.  These  areas  include  chemical  reactions  and 
particle  impact  as  they  affect  ablation.  The  experimental  effort 
is  conveniently  divided  into  two  separate  programs: 

(1)  The  chemical  erosion  program;  the  investigation  of 
chemical  effects  in  the  ablation  process  and  the  development  of 
a  simple  materials  testing  technique. 

(2)  The  particle  impact  program;  the  investigation  of  the 
effects  of  particle  impingement  on  ablation. 

As  indicated  above,  a  secondary  purpose  of  the  chemical 
erosion  program  is  to  investigate  a  relatively  simple  materials 
testing  technique.  The  usual  materials  testing  techniques  have 
taken  three  general  forms:  full-scale  rocket  tests,  small-scale 
rocket  tests,  and  splash  tests  on  materials  samples.  The  last 
technique  must  be  regarded  as  only  qualitative  since  only  a  rough 
comparison  between  the  performance  of  different  materials  is 
possible.  With  the  advent  of  large  rocket  motors,  full-scale 
rocket  tests  for  materials  evaluation  become  prohibitive.  Small- 

scale  rocket  tests  do  not  permit  the  flexibility  of  varying  test 
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parameters  independently  and  are  also  relatively  expensive.  A 
flexible,  inexpensive  yet  quantitative  experimental  materials 
evaluation  technique  is  therefore  desirable,  and  its  development 
is  the  subject  of  one  part  of  the  experimental  effort. 

The  presentation  which  follows  covers  the  areas  of  experi¬ 
mental  activity  during  the  second  quarter  of  the  project.  The 
properties  of  the  gas  mixtures  being  used  in  both  programs  were 
calculated,  and  the  results  are  given  in  tabular  and  graphical 
form.  Progress  on  the  chemical  erosion  program  is  disucssed  in 
detail;  this  includes  plasma  generator  facility  modification, 
test  section  fabrication,  and  several  firings.  The  particle 
impact  program  is  still  in  the  planning  stage,  and  the  status 
of  these  plans  is  presented. 

3.2  Properties  of  Test  Gas  Mixtures 

The  chemical  erosion  test  program,  the  experimental  investi¬ 
gation  of  specific  chemical  reactions  as  they  affect  the  erosion 
of  nozzle  wall  materials,  requires  a  series  of  gas  mixtures  that 
will  "isolate"  the  chemical  reactions  of  interest  and,  for  direct 
application  to  rocket  technology,  that  will  also  simulate  the 
solid-propellant  rocket  environment.  The  requirements  for  satis¬ 
fying  these  conditions  were  developed  and  discussed  in  Reference  1. 
In  summary,  these  requirements  were  the  duplication  of  the 
temperature-enthalpy  variation  of  the  rocket  exhaust  gases,  dupli¬ 
cation  of  the  mass  fractions  of  the  chemically  reactive  species  of 
interest,  and  duplication  of  the  heat,  mass,  and  momentum  fluxes 
at  the  surface  of  the  ablating  material. 

Five  gas  mixtures  were  selected  for  the  investigation  and  are 
summarized  below: 

(1)  The  actual  combustion  products  (with  no  aluminum  com¬ 
pounds)  . 

(2)  The  combustion  products  composition  except  for  hydrogen 
chloride  which  is  replaced  by  an  inert  gas. 
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(3)  The  mass  fractions  of  hydrogen  and  oxygen  in  the  actual 
products  that  are  capable  of  reacting  with  nozzle  wall  materials 
and  appropriate  amounts  of  inert  gases  (the  hydrogen  and  oxygen 
tied  up  in  HCl  and  CO,  respectively,  are  essentially  inert). 

(4)  The  mass  fraction  of  oxygen  that  is  capable  of  reacting 
with  nozzle  wall  materials  and  appropriate  amounts  of  inert  gases 
(the  oxygen  tied  up  in  CO  is  essentially  inert) . 

(5)  All  inert  gases. 

In  order  to  set  test  conditions  and  evaluate  test  data,  the 
thermodynamic  properties  and  composition  of  the  test  gas  mixtures 
must  be  known  for  the  complete  range  of  test  conditions.  Estimates 
of  the  elemental  composition  and  thermodynamic  properties  were 
presented  in  Reference  1.  During  this  report  period,  detailed 
calculations  of  thermodynamic  properties  and  molecular  composition 
were  made  for  each  gas  mixture  from  computer  programs  developed 
for  the  Vidya  IBM  1620  computer.  Calculations  were  performed  for 
a  complete  net  of  points  to  allow  construction  of  thermodynamic 
property  charts.  The  data  calculated  at  each  temperature  and 
pressure  were  enthalpy,  entropy,  frozen  specific  heat,  equilibrium 
specific  heat,  isentropic  exponent,  molecular  weight,  mole  fractions, 
and  mass  fractions. 

The  elemental  and  molecular  compositions  at  selected  tempera¬ 
tures  and  pressures  are  presented  in  Table  III  for  each  gas  mixture 
in  terms  of  the  mole  and  mass  fractions.  In  Mixtures  2  and  3,  the 
original  estimates  of  elemental  composition  (Ref.  1)  were  changed 
somewhat  to  yield  better  agreement  in  the  temperature-enthalpy 
variation  and  to  reduce  the  experimental  complexity  of  operating 
with  these  mixtures.  In  Mixture  2,  nitrogen  was  used  as  the  inert 
instead  of  argon  as  originally  anticipated.  In  Mixture  3,  helium 
was  eliminated,  leaving  only  nitrogen  as  the  inert-  Each  mixture 
and  its  relation  to  the  others  is  discussed  in  a  following  section. 
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The  temperature-enthalpy  variation  for  each  gas  mixture  is 
presented  in  Figures  16(a)  through  (e) .  The  variation  for  Mix¬ 
ture  1  is  included  in  each  figure  for  comparison;  this  mixture 
is  the  actual  products  of  combustion  for  a  typical  solid  propel¬ 
lant  (with  no  aluminum  compounds) .  In  each  case  the  enthalpy  base 
was  arbitrarily  assigned  as  H  =  0  at  2500^  K  and  20  atmospheres 
to  facilitate  comparison.  As  shown  in  the  figures,  the  agreement 
of  Mixtures  2  and  3  with  Mixture  1  is  excellent.  The  agreement 
is  good  for  Mixture  4  which  departs  somewhat  at  both  high  and 
low  temperatures.  Mixture  5  disagrees  substantially  at  higher 
temperatures,  above  about  3200°  K.  This  is  due  to  the  essen¬ 
tially  constant  specific  heat  caused  by  the  absence  of  chemical 
reaction.  The  only  reaction  possible,  the  dissociation  of  nitro¬ 
gen,  does  not  occur  to  any  significant  degree  in  the  temperature 
range  of  interest. 

For  convenience  in  data  reduction  and  in  setting  arc-plasma 
generator  operating  conditions,  thermodynamic  property  charts 
were  prepared  for  each  gas  mixture.  Master  charts  covering  the 
entire  range  of  temperature  and  pressure  were  made;  these  are 
presented  in  Figures  17  (a)  through  (e)  .  Plots  of  molecular 
weight  and  isentropic  exponent  and  charts  magnifying  the  primary 
area  of  interest  on  the  master  charts  were  also  made  for  each 
gas  mixture.  An  example  of  these  plots,  for  Mixture  3,  is  pre¬ 
sented  in  Figures  18(a)  through  (c)  (the  master  chart  for  Mix¬ 
ture  3  is  Fig.  17(c)). 

Transport  property  data  are  also  required  for  the  calcula¬ 
tion  of  heat-transfer  parameters.  Estimates  have  been  made  for 
the  particular  gas  mixtures  of  interest  from  the  available  data 
on  the  individual  gases  which  make  up  the  mixtures  (e.g.,  Ref.  19) 
A  computer  program  capable  of  calculating  transport  properties 
for  multicomponent  gas  mixtures  will  also  be  available  shortly. 
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The  utilization  of  and  information  obtainable  from  the  five 
gas  mixtures  will  be  discussed  in  the  following  section,  Chemical 
Erosion  Test  Program.  One  or  more  of  these  gas  mixtures  will  also 
be  utilized  in  the  Particle  Impact  Program  which  is  discussed  in 
Section  3.4. 

3.3  Chemical  Erosion  Test  Program 
3.3.1  Introduction 

The  primary  objective  of  the  chemical  erosion  test  program 
is  the  determination  of  the  effects  of  specific  chemical  reactions 
on  the  erosion  of  rocket  materials.  The  results  obtained,  in  addi¬ 
tion  to  being  of  general  importance  will  provide  necessary  guidance 
to  the  theoretical  phase  of  the  subject  contract.  A  secondary 
objective  is  the  definition  of  an  inexpensive  technique  for  the 
quantitative,  as  opposed  to  qualitative,  evaluation  of  candidate 
materials  for  rocket  nozzles.  The  technique  is  expected  to  be  a 
powerful  tool  for  the  material  fabricator  and  the  nozzle  designer. 

The  program,  in  brief,  consists  of  an  appropriate  series  of 
tests  on  ablative  materials  utilizing  the  five  gas  mixtures  dis¬ 
cussed  in  the  previous  section.  The  Vidya  arc  plasma  generator 
supplies  the  energy  input  to  the  gas  mixtures  for  simulation  of 
the  rocket  environment.  The  test  materials  are  graphite,  silica 
phenolic,  and  graphite  phenolic.  An  axisymmetric  nozzle  configu¬ 
ration  is  used  for  all  tests.  The  heat  transfer  and  pressure 
distribution  for  a  nonablating  wall  of  the  same  geometry  as  the 
test  nozzle  is  determined  with  sets  of  heat-transfer  calibration 
nozzles  and  pressure  distribution  calibration  nozzles. 

In  order  to  increase  the  effectiveness  of  the  chemical 
erosion  program,  the  originally  proposed  program  has  been  extended 
somewhat.  This  was  a  result  of  the  analytical  work  conducted 
early  in  the  program  in  which  it  was  found  that  wall  recession 
rates  are  affected  in  a  major  way  by  several  chemical  reactions. 

The  number  of  test  gases  was  therefore  increased  from  three  to 
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five,  the  test  configuration  was  changed  from  mostly  channel  tests 
to  all  nozzle  tests,  and  the  total  number  of  tests  was  increased 
from  36  to  56.  Channel  tests  were  eliminated  because  of  the  edge 
effects  that  occur  at  the  side  walls  of  the  channel  and  because  of 
the  "gouging"  which  occurs  in  the  test  material  at  the  junction 
between  the  test  material  section  and  the  water-cooled  inlet  tran¬ 
sition  section.  Also,  the  nozzle  configuration  is  directly  com¬ 
patible  with  the  axisymmetric  computer  program  being  developed. 

The  activities  during  the  report  period  that  were  directly 
related  to  the  chemical  erosion  program  included  the  following: 

(1)  Calculation  of  the  thermodynamic  properties  and  compo¬ 
sition  of  the  test  gas  mixtures. 

(2)  Completion  of  the  necessary  test  facility  changes  and 
additions . 

(3)  Checkout  of  arc-plasma  generator  operation  with  the 
test  gas  mixtures,  particularly  Mixtures  4  and  5. 

(4)  Fabrication  of  the  pressure  distribution  calibration 
nozzles . 

(5)  Redesign  and  fabrication  of  the  heat-transfer  cali¬ 
bration  nozzles. 

(6)  Procurement  of  graphite  and  silica  phenolic  nozzle 
test  section  materials. 

(7)  Fabrication  of  six  graphite  nozzle  test  sections. 

(8)  Determination  of  thermocouple  instrumentation  tech¬ 
nique  for  the  test  nozzles. 

(9)  Performance  of  a  pressure  distribution  calibration 
test  series. 

(10)  Performance  of  a  graphite  nozzle  firing. 

The  Vidya  sponsored  calibration  of  the  arc-plasma  generator 
was  also  completed  during  the  report  period.  These  activities  are 
discussed  in  following  sections. 

The  originally  anticipated  testing  schedule  for  the  report 
period  was  not  met  because  of  the  delay  in  completing  the  facility 
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calibration,  the  time  required  to  effect  the  necessary  changes 
and  additions  to  the  facility,  and  problems  in  achieving  the 
desired  arc-plasma  generator  operating  conditions.  These  diffi¬ 
culties,  discussed  in  a  following  section,  were  resolved  late  in 
the  period  and  testing  was  begun.  Pressure  distribution  tests 
were  run  with  Mixtures  4  and  5  and  a  graphite  nozzle  firing  was 
performed  with  Mixture  4.  The  results  of  these  tests  are  pre¬ 
sented  in  a  following  section.  The  chemical-erosion  program 
testing  will  be  completed  during  the  next  quarter. 

3.3.2  Program  outline 

The  program  test  matrix  is  presented  in  Table  IV.  The 
test  variables  are  the  test  gas  mixture,  the  test  section  material, 
and  the  operating  (plenum)  conditions.  The  variation  of  plenum 
pressure  is  accomplished  by  changing  the  nozzle  throat  diameter. 

A  complete  set  of  tests  is  not  run  with  Mixture  1  since  little 
or  no  difference  is  anticipated  between  it  and  Mixture  2  and 
because  of  the  hydrogen  chloride  in  Mixture  1  presents  a  hazard 
to  both  equipment  and  personnel.  Three  total  temperatures  are 
considered  since  chemical  reaction  rates  are  strongly  dependent 
on  temperature.  The  emphasis  is  on  3500°  K,  however,  since  this 
is  the  approximate  total  temperature  for  burning  of  many  solid 
propellants.  Since  reaction  rates  are  somewhat  pressure-dependent, 
four  tests  at  moderate  pressure,  175  psia,  are  being  run.  The 
majority  of  tests  are  being  run  at  300  psia  initial  plenum  pres¬ 
sure,  however.  This  value  was  chosen  as  the  upper  limit  for 
testing  in  order  to  avoid  an  unreasonably  small  nozzle  throat 
diameter  and  because  high  pressure  aggravates  the  electrode 
oxidation  problem.  The  test  section  throat  diameter  for  the 
300  psia  test  condition  is  0.3  inch  and  for  175  psia,  0.4  inch. 

For  each  gas  mixture  except  Mixture  1,  two  tests  are  conducted 
at  one  test  point  -  graphite,  3500c  K,  and  300  psia  -  to  check 
reproducibility.  Seven  tests  have  been  left  open  to  allow 
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addition  of  conditions  found  desirable  from  the  other  tests  and 
to  allow  rechecks  of  other  results  as  necessary. 

The  chemical  reaction  information  obtainable  from  the  test 
program  may  be  discussed  conveniently  by  referring  to  the  molecular 
compositions  of  Table  III.  The  influence  of  hydrogen  chloride  on 
the  wall  material  erosion  is  obtainable  from  comparison  of  results 
from  Mixtures  1  and  2.  As  mentioned  previously,  HCl  is  expected 
to  act  as  an  inert  with  the  wall  materials.  In  going  from  Mix¬ 
ture  2  to  Mixture  3,  the  carbon  compounds,  CO  and  CC>2 ,  are  elimi¬ 
nated.  The  effect  of  their  removal  on  wall  erosion  is  assessed 
here.  The  major  change  is  that  the  available  oxygen  in  the  form 
of  C02  in  Mixture  2  is  found  in  the  form  of  water  in  Mixture  3. 

The  CO  is  expected  to  act  as  an  inert  and  therefore  its  removal 
may  not  have  an  affect  on  the  wall  erosion. 

In  going  from  Mixture  3  to  Mixture  4,  hydrogen  is  elimi¬ 
nated.  The  oxygen  is  in  the  form  of  02 ,  O,  and  NO  rather  than  as 
mostly  H20  as  in  Mixture  3.  This  shift  in  the  form  of  oxygen  and 
the  possible  affect  of  atomic  hydrogen  recombination  on  the  corro¬ 
sion  of  the  wall  material  will  be  studied  here.  Mixture  5  will 
serve  to  determine  the  ablation  in  the  absence  of  chemical  reaction 
with  the  free  stream.  In  this  case,  the  mechanical  properties  and 
thermal  degradation  characteristics  of  the  material  will  deter¬ 
mine  the  surface  recession. 

3.3.3  Test  facility  modifications  and  checkout 

A  number  of  changes  and  additions  to  the  plasma  generator 
facility  were  required  to  accommodate  the  chemical  erosion  test 
program.  These  modifications  were  primarily  associated  with  the 
large  number  of  gases  being  used  and  the  safety  and  operating 
problems  associated  with  some  of  them.  Separate  gas  feed  and 
metering  systems  were  built  for  hydrogen,  carbon  dioxide,  and 
hydrogen  chloride,  which  must  be  introduced  into  the  plasma  gene¬ 
rator  separately.  Another  system  was  built  for  the  mixtures  of  the 
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other  gases  that  are  to  be  introduced  as  the  primary  plasma- 
generator  gas.  A  schematic  of  the  gas  system  is  presented  in 
Figure  19.  The  necessary  safety  precautions  were  observed  in 
building  up  each  system.  A  gas  collection  and  exhaust  system 
was  also  installed  close  to  the  exit  of  the  plasma  generator  to 
evacuate  the  sometimes  toxic  and  explosive  exhaust  gases  from 
the  laboratory. 

An  electrically  heated  water  bath  was  constructed  to  heat 
the  condensable  gases  (C02  and  HC1)  in  order  to  eliminate  severe 
gas  bottle-pressure  drops  and  condensation  problems  when  drawing 
from  these  bottles.  A  temperature  of  approximately  110°  F  is 
maintained  in  the  bath. 

Checkout  tests  were  run  on  the  arc-plasma  generator  to 
determine  the  operating  characteristics  and  conditions  for  the 
gases  of  interest.  Because  of  the  low  molecular  weights  of  the 
primary  gas  mixtures  (the  gases  which  are  arc  heated) ,  arc  insta¬ 
bility  problems  were  encountered.  The  combination  of  low  molec¬ 
ular  weight,  high  pressure,  and  the  relatively  low  enthalpy, 
required  high- operating  voltage  which  was  incompatible  with  the 
power  supply.  Instability  and  eventual  blow-out  of  the  arc 
occurred.  The  problem  was  alleviated  by  bypassing  some  of  the 
primary  gas  and  introducing  it  downstream  of  the  arc.  The 
resultant  decrease  in  gas-flow  rate  through  the  arc  and  in  arc 
voltage  resulted  in  stable  and  smooth  operation.  Further  tests 
were  run  to  determine  the  required  preset  conditions  to  achieve 
the  desired  test  conditions  for  Mixtures  4  and  5. 

The  Vidya  sponsored  calibration  of  the  arc-plasma  gene¬ 
rator  was  also  completed  during  the  report  period.  The  results 
are  presented  in  Reference  20  and  are  discussed  briefly  here. 

A  compact  heat-exchanger  calorimeter,  schematically  shown  in 
Figure  20,  was  connected  to  the  exit  of  the  plasma  generator  to 
remove  the  energy  in  the  exit  gas  stream.  The  energy  input  to 
the  arc-heated  gas  as  calculated  by  the  conventional  energy- 
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that  absorbed  by  the  calorimeter.  The  results  are  presented  in 
Figure  21;  almost  all  test  points  fall  within  a  +15  percent  band. 
The  results  of  this  program  indicate  that  test  conditions  calcu¬ 
lated  by  the  energy  balance  method  are  reasonably  accurate  and 
that  the  test  results  are  reproducible. 

3.3.4  Test  sections 

3. 3. 4.1  Calibration  nozzles 

Pressure  distribution  and  heat-transfer  calibration 
nozzles  were  fabricated  to  determine  the  nonablating  wall-pressure 
distribution  and  heat  transfer  for  a  nozzle  geometry  identical  to 
that  of  the  material  nozzle  configurations.  The  data  obtained 
from  tests  with  these  nozzles  are  required  for  the  analysis  and 
interpretation  of  the  data  obtained  in  the  material  nozzle  firings. 

As  discussed  in  Reference  1,  the  pressure  distribution 
nozzles  are  uncooled;  the  heat-sink  capabilities  of  the  copper 
bodies  allow  sufficient  run  times  to  obtain  a  steady-state  pressure 
distribution.  A  photograph  of  one  of  the  nozzles  installed  on  the 
plasma  generator  is  shown  in  Figure  22.  The  pressure-tap  locations 
are  presented  in  Figure  23.  The  pressure  distribution  is  measured 
with  a  high-pressure  manometer  and  pressure-gauge  board  in  which 
the  pressure  readings  are  "trapped"  by  simultaneously  closing  a 
series  of  quick-acting  valves  after  steady-state  is  reached.  The 
board  is  then  read  and  data  recorded  after  the  test.  The  board 
is  visible  on  the  right  in  Figure  22. 

The  set  of  calibration  nozzles  consists  of  three  units 
of  different  throat  diameter,  0.3,  0.4,  and  0.5  inch.  This  range 
is  required  to  cover  both  the  pressure  conditions  at  which  tests 
are  being  run  (see  Table  IV)  and  to  provide  information  on  the 
pressure  distribution  in  the  ablative-material  nozzles  as  the 
throat  diameter  increases  during  each  test. 

Because  of  a  fabrication  problem,  the  original  design 
of  the  heat-transfer  calibration  nozzles  (Ref.  1)  was  changed. 
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Brief  ly,  the  fabrication  technique  was  the  electroforming  of  a 
copper-nozzle  insert  which  contained  thermocouple  instrumentation 
for  determination  of  the  heat-flux  and  surface-wall  temperature. 

The  thermocouples  were  set  in  position  before  electroforming  was 
begun  and  were  therefore  an  integral  part  of  the  nozzle  insert. 
However,  no  thermocouple  insulation  or  protective  coating  was 
found  that  would  holdup  under  the  many  hours  of  immersion  in  the 
electrolyte  during  the  plating  process.  Also,  due  to  the  differ¬ 
ence  in  the  thermal  conductivity  of  constantan  and  copper,  heat 
conduction  normal  to  the  surface  could  cause  significant  errors 
in  the  measured  surface  temperatures.  The  original  design  was 
therefore  discarded  in  favor  of  a  steady-state  calorimetric  mea¬ 
surement  of  heat  transfer,  a  method  which  has  been  used  success¬ 
fully  in  the  past  (e.g.  Refs.  21  and  22). 

The  nozzle  configuration  is  presented  in  Figure  24. 

It  is  madeup  of  individually  water-cooled  segments  which  are 
stacked  together  to  form  the  nozzle.  The  throat  region  is,  of 
course,  the  region  of  most  interest  and  therefore  contains  the 
most  segments.  The  nozzle  wall-material  is  copper.  A  0.002-inch 
gap  between  segments  (Fig.  24)  serves  as  an  insulator  between  the 
heat-transfer  surfaces  of  adjacent  segments.  The  water-flow 
rates  in  each  segment  are  monitored  by  rotameters  and  the  thermo¬ 
couple  outputs  are  recorded  on  an  oscillograph.  The  set  of  three 
calorimeter  nozzles  have  throat  diameters  of  0.3,  0.4,  and  0.5  inch, 
the  same  as  the  pressure-distribution  nozzles. 

The  heat-flux  for  each  segment  is  determined  from  its 
measured  water-flow  rate  and  temperature  rise.  With  the  heat 
input  to  the  cooling  water  calculated  from  these  measurements 
and  the  known  heated-surface  area  for  each  segment,  the  local 
average  heat-flux  is  then  calculated.  This  heat-flux  is  due  to 
both  convective  and  radiative  heating.  The  convective  heat-flux 
is  expected  to  be  the  greater  of  the  two  in  all  regions  of  the 
nozzle.  The  radiative-flux  is  due  to  radiation  interchange 
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between  the  gas  stream  and  the  nozzle  wall  and  between  the  arc 
and  the  wall.  This  contribution  to  the  total  heat-flux  is 
expected  to  be  small  in  the  absence  of  particles.  The  experi¬ 
mental  results  will  be  compared  with  theoretical  predictions 
for  convective  heating  (i.e.,  Bartz ,  ARM).  Estimates  of  the 
radiative-flux  will  also  be  made. 

3. 3. 4. 2  Ablative-material  nozzles 

The  nozzle  firing  phase  of  the  program  was  discussed 
in  some  detail  in  Reference  1.  In  brief,  the  test  section 
materials  are  graphite,  silica  phenolic,  and  graphite  phenolic. 

The  configuration  of  the  0.3-inch  throat  diameter  nozzle  is 
presented  in  Figure  25  for  convenient  reference.  The  internal 
wall-temperature  distribution  is  measured  as  a  function  of  time 
to  allow  determination  of  hot-side  wall  temperature.  The  deter¬ 
mination  of  surface  recession  is  made  throughout  the  entire 
nozzle  from  pre-  and  post-test  measurement  of  the  internal-wall 
profile  and  at  the  throat,  (the  downstream  end  of  the  region  of 
minimum  cross-sectional  area,  as  discussed  in  Section  3. 3. 5.1), 
from  calculations  based  on  the  measured  time  variation  of  plenum 
pressure . 

Because  of  its  widespread  usage,  ATJ  graphite  was  chosen 
as  the  graphite-nozzle  material.  Some  of  its  shortcomings  for 
this  application  have  become  apparent  in  the  meantime  however. 

After  test,  the  wall  surface  in  the  throat  region  of  the  first 
graphite  nozzle  was  quite  rough  and  pitted.  The  low  density  and 
therefore  high  porosity  of  ATJ  graphite  is  felt  to  be  at  least 
the  partial  cause  of  this  surface  condition.  The  condition  could 
be  the  result  of  spalling  or  may  be  due  to  nonuniformities  in 
porosity  which  cause  nonuniform  material  removal  by  mechanical 
and/or  chemical  action.  The  high  porosity  could  allow  the  material 
to  absorbe  significant  amounts  of  moisture  or  other  gases.  When 
the  material  is  heated  at  high  heat-flux,  these  gases  may  cause 
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spalling  due  to  the  rapid  expansion  of  the  gases.  A  further  con¬ 
sequence  of  the  high  porosity  is  that  the  surface  condition, 
microscopically,  is  quite  rough  and  therefore  offers  a  larger 
exposed  surface  area  than  the  actual  planer  surface  and  a  large 
number  of  active  sites  for  chemical  reaction. 

Because  of  these  problems ,  other  higher-density  graphites 
were  investigated.  Graph-i-tite ,  grade  GX,  was  determined  to  be 
a  reasonable  possibility;  it  has  a  high  density  (Sp.  gr.  1.90  as 
compared  to  1.73  for  ATJ) ,  good  thermal-shock  characteristics, 
and  is  reasonably  isotropic.  A  length  of  this  graphite  was  there¬ 
fore  procured  and  a  preliminary  test  will  be  run  to  check  its 
performance.  A  decision  will  then  be  made  on  the  graphite  type 
to  be  used  for  all  following  tests. 

The  internal-temperature  distributions  of  the  test  nozzles 
are  being  determined  in  order  to  allow  calculation  of  the  wall- 
surface  temperature  desired  for  the  interpretation  of  test  results. 
(Data  reduction  will  be  accomplished  with  the  internal  conduction 
program.)  Originally,  a  plug  technique  was  chosen  as  the  method 
of  thermocouple  instrumentation  to  measure  these  distributions. 
Problems  in  fabrication  and  possible  difficulties  in  data  inter¬ 
pretation  dictated  a  change  in  this  decision.  The  machining  of 
small  tapered  or  threaded  holes  and  the  machining  of  the  plugs 
themselves  was  found  to  be  extremely  difficult  and  time  consuming 
for  the  phenolic  test  sections.  The  data-interpretation  problem 
was  due  to  the  unknown  effect  of  contact  resistance  at  the  inter¬ 
face  between  the  plug  and  nozzle  body.  The  plug  configuration 
does  not  coincide  with  the  material  radial- direction  of  heat-flux 
and  therefore  any  contact  resistance  at  the  nozzle  body  -  plug 
interface,  results  in  a  distortion  of  the  heat  flux  (see  following 
sketch).  A  reliable  estimate  of  this  contact  resistance  for  use 
in  correcting  for  this  effect  is  difficult,  if  not  impossible.  Hence, 
other  thermocouple  instrumentation  techniques  were  investigated. 


A  prime  requirement  in  thermocouple  instrumentation,  of 
course,  is  that  the  sensor  does  not  disturb  the  heat  flow  at  the 
point  of  measurement.  To  satisfy  this  requirement,  the  sensor 
should  be  as  small  as  practical  and  the  thermocouple  leads  in  the 
vicinity  of  the  junction  should  be  placed  along  an  isotherm.  This 
last  consideration  is  particularly  important  in  low  conductivity 
materials  and  is  impressively  illustrated  in  Reference  23.  Also, 
the  thermocouple  junction  must  be  in  intimate  contact  with  the 
material  in  order  to  eliminate  errors  due  to  contact  resistance. 
Another  obvious  requirement  is  that  the  location  of  the  junction 
be  precisely  known.  After  considering  these  requirements  and  the 
practical  aspects  of  fabrication,  the  thermocouple  installation 
scheme  shown  in  Figure  26  was  selected.  A  photograph  of  one  of 
the  thermocouple  probes  is  presented  in  Figure  27.  The  probes 
are  inserted  into  the  test  nozzle  so  that  they  are  tangent  to  an 
isotherm  at  their  junction.  The  diameter  of  the  ceramic  insulator 
which  encloses  the  thermocouple  wires  is  0.035  inch  and  therefore 
offers  a  minimum  disturbance  to  the  heat  flux.  The  thermocouples 
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are  spring  loaded  against  the  bottom  of  the  thermocouple  probe 
hole  to  insure  good  contact  of  the  thermocouple  junction  with 
the  wall.  The  thermocouple  wires  are  36  gauge  (0.005-inch-diameter). 
Chrome 1-Alumel  thermocouples  are  used  where  temperatures  are  not 
expected  to  exceed  2500°  F.  Thermocouples  of  tungsten  5  percent 
rhenium  and  tungsten  26  percent  rhenium  are  used  for  higher  tempera¬ 
ture  locations.  The  small  gauge  of  the  wires  further  reduces  the 
possibility  of  error  due  to  conduction  along  the  wires.  The  loca¬ 
tions  of  the  thermocouple  junctions  are  determined  by  x-ray  photog¬ 
raphy.  The  requirements  for  obtaining  accurate  data,  discussed 
above,  are  therefore  closely  satisfied. 

The  thermocouple  instrumentation  is  located  at  three  axial 
positions  in  the  tubular  throat  as  shown  in  Figure  25.  A  typical 
installation  at  one  of  these  positions  is  illustrated  in  Figure  26. 
The  thermocouple  radial  locations  and  the  analytically  predicted 
temperature  distributions  for  test  nozzles  of  each  material  are 
presented  in  Figures  28(a)  -  (c) .  The  distributions  are  idealized 
in  that  room  temperature  thermal  properties  were  used,  no  ablation 
was  assumed,  and,  for  the  phenolics,  no  allowance  was  made  for  char 
formation.  The  temperature  distributions  are  for  the  throat  region 
of  a  0.3-inch  diameter  nozzle  and  were  obtained  from  the  tabulated 
solutions  of  Reference  24  for  convective  heating  of  a  hollow 
cylinder.  The  results  are  for  the  highest  heat-flux  test  condition 
(Tq  =  3500°  K,  pQ  =  300  psia) .  The  thermocouple  locations  were 
chosen  from  consideration  of  these  temperature  distributions  and 
the  expected  ablation  rates. 

3.3.5  Test  results 

Testing  during  the  report  period  that  was  a  direct  part  of 
the  chemical  erosion  program  consisted  of  a  series  of  pressure 
distribution  calibration  tests  and  an  ATJ  graphite  nozzle  firing. 

The  pressure  distribution  tests  were  performed  with  Mixtures  4  and 
5  at  the  high-temperature  test  conditions  and  with  the  0.3-inch 
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and  0.4-inch  throat  diameter  nozzles.  These  conditions  correspond 
to  Test  Points  28-31,  32,  39-42,  and  43  (see  Table  4).  The  nozzle 
firing  was  performed  at  Test  Point  28  -  Mixture  4  at  the  nominal 
conditions  3500°  K  total  temperature  and  300  psia  initial  plenum 
pressure.  No  thermocouple  instrumentation  was  included  in  this 
test  because  of  its  preliminary  nature  as  the  first  graphite 
nozzle  firing. 

3. 3. 5.1  Pressure  distribution  tests 

The  pressure  distribution  results  are  presented  in  Fig¬ 
ures  29(a)  and  (b)  and  are  discussed  below.  The  test  conditions 
for  each  test  are  indicated  on  the  figures.  In  the  tests  with 
Mixture  5,  the  enthalpy  and  therefore  temperature  were  inadvertently 
run  higher  than  planned. 

The  circumferential  variation  of  pressure  was  measured  at 
the  central  section  of  the  throat  with  four  pressure  taps  90°  apart. 
In  all  tests,  the  variation  was  less  than  2  psi  which  corresponds 
to  a  variation  of  less  than  1  percent  of  the  stagnation  pressure. 

A  fifth  test,  not  presented  in  the  figures,  was  run  with  Mixture  4 
on  the  0.4-inch  diameter  throat  nozzle  to  check  reproducibility. 

The  points  for  this  test  fell  almost  identically  on  the  appropriate 
points  of  Figure  29(b). 

The  theoretical  pressure  distribution  for  an  idealized 
one-dimensional  isentropic  expansion  in  the  nozzle  are  also  pre¬ 
sented  in  Figures  29.  In  a  constant  diameter  tubular  throat,  the 
theoretical  pressure  in  the  absence  of  friction  and  heat  transfer 
is  constant.  Actually  however,  boundary-layer  effects  will  cause 
a  pressure  drop.  The  boundary  layer  would  be  expected  to  grow 
with  increasing  distance  down  the  throat  and  therefore  a  decreasing 
effective  flow  area  would  result.  The  sonic  point,  the  point  of 
minimum  area,  would  then  be  at  the  throat  exit.  If  the  tubular 
throat  diameter  were  not  exactly  constant  throughout  due  to  machin¬ 
ing  nonuniformities,  the  pressure  drop  through  the  throat  and  the 
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location  of  the  sonic  point  would  be  affected  further.  The  combined 
effects  of  these  nonuniformities  and  boundary- layer  growth  would 
result  in  an  appreciable  pressure  drop  and  the  sonic  point  would 
be  expected  somewhat  downstream  of  the  minimum  physical  diameter 
in  the  throat  region. 

Returning  now  to  Figures  29,  the  agreement  of  the  experi¬ 
mental  results  with  the  predicted  distribution  is  good  in  the  con¬ 
traction  and  expansion  sections.  In  the  tubular  throat,  however, 
the  agreement  is  not  good  and,  based  on  the  above  discussion,  is 
not  too  surprising.  A  constant  pressure  is  predicted  by  theory 
whereas  a  pressure  drop  of  approximately  58  psi  occurs  for  the 
0.3-inch  diameter  throat  and  approximately  18  psi  for  the  0.4-inch 
throat.  This  disagreement  can  be  attributed  to  machining  nonuni¬ 
formities  and  boundary- layer  effects  discussed  above.  Variations 
occur  in  the  real  and  effective  flow  area  through  the  tubular 
throat  which  have  a  significant  effect  on  the  pressure  distri¬ 
bution. 

Measurements  of  the  variation  in  diameter  between  the 
pressure  taps  in  the  throat  (2  and  7)  were  made  by  traversing  the 
throat  with  a  pivot-arm  dial  gauge  micrometer.  For  the  0.3-inch 
throat,  the  variation  was  0.0046-inch  over  the  distance  between 
taps  and  the  minimum  diameter  occurred  close  to  the  entrance  to 
the  throat.  The  pressure  drop  corresponding  to  the  resultant 
variation  in  area  is  approximately  39  psi.  The  boundary- layer 
effect  was  also  investigated  by  a  simple  calculation  of  the  pres¬ 
sure  drop  due  to  the  combined  effects  of  friction  and  heat  transfer 
(Ref.  25).  For  the  0.3-inch  throat,  this  contribution  is  approxi¬ 
mately  13  psi.  The  total  predicted  pressure  drop  due  to  machining 
nonuniformities  and  boundary-layer  effects  is  therefore  52  psi  as 
compared  to  the  measured  value  of  58  psi.  The  agreement  is  sur¬ 
prisingly  good.  As  seen  in  Figure  29(a) ,  the  experimental  sonic 
points  for  both  tests  are  approximately  2/3  of  the  distance  down 
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the  tubular  throat.  The  combined  effects  of  machining  nonuni¬ 
formities  and  boundary- layer  growth  apparently  resulted  in  a  mini¬ 
mum  effective  area  at  this  point. 

A  similar  analysis  was  made  for  the  0.4-inch  throat 
nozzle.  The  variation  in  throat  diameter  was  0.0020-inch  and 
results  in  a  pressure  drop  of  approximately  13  psi.  The  combined 
effects  of  friction  and  heat  transfer  contribute  an  8  psi  drop. 

The  total  predicted  drop  here  is  then  21  psi  as  compared  to  the 
measure  value  of  18  psi.  Again  the  agreement  is  surprisingly  good. 
The  minimum  physical  throat  diameter  occurred  approximately  half¬ 
way  down  the  throat.  The  sonic  point,  as  seen  from  Figure  29(b), 
is  close  to  the  throat  exit.  Again  boundary- layer  growth  apparently 
caused  the  effective  area  to  be  a  minimum  downstream  of  the  point 
of  the  physical  minimum  area. 

The  experimental  results  indicate  that  the  pressure 
distributions  in  the  contraction  and  expansion  sections  of  the 
nozzles  are  closely  defined  by  a  one-dimensional  isentropic  ex¬ 
pansion.  This  does  not  apply  in  the  tubular  throat,  but  the 
pressure  drop  here  can  be  predicted  quite  closely  by  consider¬ 
ing  both  the  actual  nonuniformities  in  throat  diameter  and 
boundary- layer  effects.  Further  pressure  distribution  tests  will 
be  performed  to  fully  confirm  this;  testing  over  the  complete 
range  of  test  conditions  is  not  anticipated  because  of  the  favor¬ 
able  results  obtained  so  far. 

3. 3. 5. 2  Chemical  erosion  test 

The  ATJ  graphite-nozzle  firing  results  are  presented  in 
Figures  30  and  31.  A  photograph  of  the  firing  is  shown  in  Fig¬ 
ure  32  and  a  post-test  photograph  of  the  nozzle  is  presented  in 
Figure  33. 

Figure  30  shows  the  before  and  after  internal  profile 
of  the  nozzle  at  one  cross-section.  The  maximum  surface  reces¬ 
sion  occurred  close  to  the  entrance  to  the  throat  and  corresponded 
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to  an  average  recession  rate  over  the  31.0  seconds  of  test  of 
approximately  6.2  mils/sec.  Over  most  of  the  tubular  throat 
region,  the  recession  was  approximately  uniform  at  an  average 
rate  of  about  5.2  mils/sec.  The  recession  was  sufficiently  great 
that  no  expansion  section  was  left  at  the  end  of  test;  the  minimum 
diameter  was  at  the  exit  of  the  test  section.  The  circumferential 
uniformity  of  erosion  was  difficult  to  measure  in  the  throat  region 
because  of  the  roughness  of  the  final  surface,  however  visual 
inspection  indicated  no  significant  assymetry.  Measurements  were 
made  in  the  contraction  cone  and  the  uniformity  was  good.  The  vari¬ 
ation  at  one  plane  between  two  mutually  perpendicular  measurements 
of  diameter  did  not  exceed  0.008  inch  or  1.5  percent  of  the  dia¬ 
meter  . 

After  test,  the  entire  surface  in  the  tubular  throat 
region  was  quite  rough  as  seen  in  Figure  33.  As  mentioned  pre¬ 
viously,  this  condition  is  felt  to  be  at  least  partially  due  to 
the  low  density  and  therefore  high  porosity  of  ATJ  graphite. 

Because  of  this  condition,  it  appears  that  the  mechanical  removal 
of  material  was  a  significant  contributor  to  the  surface  recession. 

A  higher  density  graphite  will  be  tested  shortly  and  a  choice  will 
then  be  made  as  to  the  graphite  type  to  be  used  in  the  rest  of  the 
test  program. 

The  instantaneous  throat  surface  recession  may  be  deter¬ 
mined  for  any  time  during  a  test  firing  from  the  time  history  of 
plenum  pressure.  The  instantaneous  nozzle  throat  area  is  directly 
related  to  the  plenum  pressure  through  the  equation 
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where  C  is  a  nozzle  coefficient  (C  <  1)  which  accounts  for 
non-isentropic  effects  due  to  friction  and  heat  transfer.  The 
equation  assumes  that  both  the  molecular  weight  M  and  isentropic 


exponent  y  are  constant  in  the  expansion  process  from  plenum  to 
sonic  conditions.  For  the  test  gas  mixtures,  these  properties 
actually  vary  somewhat  in  the  expansion.  Their  variation  is  small, 
however,  and  therefore  the  equation  yields  accurate  results. 

Sufficient  data  are  taken  in  each  nozzle  firing  to  deter¬ 
mine  all  variables  of  Equation  (63) .  The  nozzle  coefficient  C 

n 

is  determined  as  that  value  which  makes  the  measured  throat  area 
before  test  equal  to  the  calculated  throat  area  at  the  start  of 
test  before  surface  recession  has  begun.  The  throat  area  and  dia¬ 
meter  (A*  and  D*)  are  then  calculated  as  a  function  of  time  from 
the  measured  time  variation  of  plenum  pressure  and  other  operating 
conditions.  The  surface  recession,  assuming  the  location  of  the 
throat  does  not  shift  with  time,  is  then  given  at  any  time  0 
simply  by 


0  = 


<D*>e 


(D*>e 


(64) 


The  results  of  this  calculation  for  the  graphite  nozzle 
firing  are  presented  in  Figures  31(a)  through  (c) ;  the  plenum  pres¬ 
sure,  the  surface  recession,  and  the  surface  recession  rate  are 
presented  as  functions  of  test  time.  Unfortunately,  as  seen  in 
Figure  30,  the  nozzle  throat  moved  downstream  with  increasing  test 
time  until,  at  the  end  of  test,  the  throat  was  actually  at  the 
nozzle  exit.  The  evaluation  of  surface  recession  based  on  the 
throat  calculations  are  therefore  lower  than  the  actual  recession 
in  what  was  the  throat  region  before  test.  To  eliminate  this  prob¬ 
lem  in  future  tests,  consideration  is  being  given  to  eliminating 
the  expansion  section  in  the  test  nozzles  and  decreasing  the  heat¬ 
sink  effect  of  the  nozzle  retaining  ring. 

Turning  now  to  the  test  results,  Figure  31(b)  shows  the 
surface  recession  as  a  function  of  test  time.  It  can  be  seen  that 
the  transient  time  to  reach  a  quasi-steady  ablation  condition  is 
approximately  8  seconds .  The  total  predicted  surface  recession 
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(using  Eq.  (63))  after  the  31.0  seconds  of  test  is  102  mils.  This 
corresponds  to  a  final  throat  diameter  of  0.513  inch  as  compared 
to  the  measured  value  after  test  of  0.554  inch.  This  discrepancy 
may  well  be  due  to  the  rough  surface  condition  causing  a  decrease 
in  the  nozzle  coefficient,  Cn ,  during  the  firing.  Also,  the 
minimum  diameter,  the  sonic  point  for  which  Equation  (63)  is  valid, 
moved  downstream  with  increasing  test  time  until,  at  the  end  of 
test,  it  was  at  the  nozzle  exit.  The  calculated  recession  is  there¬ 
fore  less  than  that  for  the  tubular  throat  region.  As  discussed 
previously,  efforts  are  being  made  to  alleviate  this  problem. 

The  instantaneous  surface-recession  rate  determined  from 
Figure  31(b)  is  presented  in  Figure  31(c).  The  trend,  a  decreasing 
recession  rate  with  increasing  time,  is  as  expected.  An  increase 
in  diameter  results  in  a  decrease  in  the  heat  flux  and  the  rate  of 
species  diffusion  to  the  wall,  thus  decreasing  the  chemical  erosion 
rate . 

The  theoretical  surface-recession  rate  was  also  calculated 
from  the  dimensionless  wall-recession  rate  for  Mixture  4  predicted 
by  Equation  (53) . 


This  equation  assumes  that  chemical  equilibrium  exists  at  the  wall 
and  that  chemical  reaction  is  the  only  wall  erosion  mechanism. 

Recall  that  Equation  (53)  is  essentially  independent  of  wall  tempera¬ 
ture  for  Mixture  4  in  the  temperature  range  2000'"'  K  to  3300c  K  (see 
Section  2.3.4),  thus  obviating  the  need  for  coupling  the  surface 
chemistry  considerations  with  the  internal  conduction  program.  From 
Table  III,  the  elemental  mass  fraction  for  Mixture  4  is  =  0.1525 

and  therefore  Equation  (53)  reduces  to 


B’  =  0.1144 


(65) 


B1  is  defined  in  Equation  (48)  as 


re  e  M 

and  is  the  ratio  of  the  rate  of  introduction  of  wall  material  into 
the  boundary  layer  to  the  rate  of  species  diffusion  to  the  wall. 
The  term  (pv)w  is  related  to  the  material  behavior  by 

(Pv)w  =  Pm  If  (6e 


where  p^  is  the  density  of  the  material  and  da/d0  is  the 
surface  recession  rate.  From  Equations  (48),  (65),  and  (66),  the 

theoretical  surface-recession  rate  for  Mixture  4  is  then  given  by 


“  °*1144 


PeueCM 


The  term  p  u  C,.  was  determined  from  the  heat-transfer  coefficient 
from  the  simplified  Bartz  Equation  (Eq.  (56),  Ref.  1).  The  Bartz 
equation  was  used  instead  of  the  ARM  method  because,  as  stated  in 
Section  2.2.1,  it  is  felt  that  the  ARM  method  requires  further 
verification  before  it  can  be  recommended  for  all  nozzle  con¬ 
ditions.  The  diameter  used  in  the  calculation  was  that  derived 
from  the  calculated  throat  area,  Equation  (63),  and  the  Lewis- 
Semenov  number  was  assumed  to  be  unity 


(p  u  C„  =  p  u  CT7  Le2/3) 
re  e  M  <=  w 


e  e  H 


The  surface-recession  rate  as  predicted  from  Equation  (67) 
is  indicated  on  Figure  31(c)  as  the  dashed  line.  The  agreement 
with  the  experimental  results  is  amazingly  good,  but  a  word  of 
caution  is  in  order.  In  comparing  results,  it  must  be  remembered 
that  the  surface  recession  presented  in  Figure  31(b)  was  lower 
than  that  measured  after  test  by  16  percent.  Based  on  this  post¬ 
test  measurement,  the  instantaneous  experimental  surface-recession 


rate  would  be  somewhat  greater  than  that  shown  by  the  solid  line 
of  Figure  31(c),  and  thus  greater  than  the  theoretical  prediction. 
This  would  be  indicative  of  mechanical  erosion  in  addition  to  the 
contribution  of  chemical  reaction  or,  as  discussed  previously,  a 
higher  effective  surface  area  for  chemical  reaction.  Also  note 
that  Equation  (67)  assumes  that  chemical  equilibrium  obtains  at 
the  wall.  If  reaction  is  kinetically  controlled,  the  contribution 
of  chemical  action  to  the  surface  recession  is  reduced,  and  the 
dashed  line  of  Figure  31(c)  would  be  lowered.  Again,  the  result¬ 
ing  difference  could  be  attributed  to  mechanical  erosion.  Further 
had  the  ARM  method  been  used  with  Equation  (67)  in  generating  the 
dashed  cruve,  the  curve  would  be  lowered  by  about  30  percent,  but 
would  retain  the  same  shape.  If  it  is  found  that  the  ARM  method 
is  valid  for  the  test  conditions  used,  the  resulting  discrepancy 
between  predicted  and  experimental  curves  would  have  to  be  attri¬ 
buted  to  mechanical  erosion. 

It  is  seen  that  the  agreement  between  experiment  and 
theory  is  not  as  clear  cut  as  might  be  thought  at  first  glance, 
primarily  because  of  the  counteracting  effects  of  mechanical 
erosion  and  reaction  kinetics.  The  test  program  is  set  up  so 
that  these  individual  contributions  may  be  assessed;  thus,  a 
complete  explanation  of  the  results  of  Figure  31(c)  will  await 
further  testing. 

3.4  Particle  Impact  Test  Program 

The  particle  impact  program,  the  other  phase  of  the  experi¬ 
mental  effort  in  the  subject  contract,  is  concerned  with  the 
effects  of  particle  impingement  on  the  material  behavior  in  a 
rocket  nozzle.  The  various  mechanisms  associated  with  the  impact 
and  presence  of  particles  are  therefore  investigated.  These 
mechanisms  include  thermal-energy  transfer,  kinetic-energy  trans¬ 
fer,  shear,  chemical  reactions,  and  effects  on  the  convective  and 
radiative  heat-transfer. 


[ 


Because  of  the  large  effort  required  by  the  chemical  erosion 
test  program,  the  progress  in  this  phase  of  the  program  has  been 
minimal.  The  primary  effort  to  date  has  been  the  specific  difi- 
nition  of  instrumentation  techniques  and  procedures.  This  instru¬ 
mentation  is,  of  necessity,  somewhat  unique  and  complicated  in 
order  to  yield  the  desired  test  data.  The  general  facility  setup 
and  checkout  in  preparation  for  the  program  has,  for  the  most  part, 
been  completed  in  that  the  similar  activities  in  preparation  for 
the  chemical  erosion  program  are  also  applicable  here.  The  program 
and  instrumentation  are  discussed  briefly  below. 

The  Vidya  arc-plasma  generator  is  utilized  to  simulate  the 
rocket-exhaust  environment.  The  test  gas  to  be  used  will  be  the 
simplest  of  the  gas  mixtures  discussed  previously  which  also 
adequately  duplicates  the  nozzle  materials  response.  Particles 
will  be  introduced  in  the  plenum  chamber  of  the  plasma  generator 
by  a  powder-hopper  feed  system  or  by  a  particle-slurry  feed  along 
with  the  secondary  injection  gases. 

The  test  configuration  consist  of  a  flat-plate  slab  of  the 
test  material  and  is  positioned  in  the  hot-gas  stream  just  down¬ 
stream  of  the  plasma  generator  nozzle  exit.  The  test  slab  mate¬ 
rials  are  tungsten,  graphite,  silica  phenolic,  and  graphite  phenolic 
Tungsten  is  included  for  two  reasons  -  first,  it  is  a  material  used 
in  rocket  motors  and  second,  it  is  a  nonablator  and,  therefore, 
serves  as  a  'clean'  basis  for  assessing  particle-impact  effects 
without  the  added  complication  of  ablation.  The  angle  of  attack 
of  the  test  slab  with  respect  to  the  exit  flow  is  variable.  Plasma 
generator  operation  will  first  be  established  either  with  or  with¬ 
out  particles  in  the  gas  stream  and  then  the  test  slab  swung  into 
position.  Particle  flow  will  be  established  before  or  at  any  time 
after  the  test  material  is  exposed  to  the  gas  stream.  Control 
tests  with  no  particle  loading  will  also  be  run  to  checkout  the 
material  performance  with  no  particles  present.  The  temperature 
response  of  the  material  will  be  monitored. 
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The  test  variables  are  the  test-slab  material,  the  slab  wall- 
temperature,  the  plasma-generator  operating  conditions  which  will 
also  influence  the  particle  temperature  and  particle  velocity,  the 
particle  size,  the  particle  mass  loading,  and  the  angle  of  inci¬ 
dence  of  the  particles  on  the  test  slab.  The  range  of  the  last 
three  parameters  will  be  dictated  in  part  by  input  from  Aeronu- 
tronic  under  a  companion  program. 

A  number  of  special  instrumentation  techniques  are  required 
for  data  acquisition.  Because  of  the  high  particle  speeds,  up  to 
5000  ft/sec,  an  image  intensifier  camera  will  be  used  for  measure¬ 
ment  of  particle  velocity  both  before  and  after  impact.  Film 
exposure  is  by  electron  bombardment  rather  than  by  the  light 
emission  from  the  particles.  The  extremely  short  exposure  times 
required  are  therefore  possible.  Velocities  will  be  determined 
by  a  streak-photography  method  similar  to  that  used  in  star  and 
satellite  tracking  in  which  an  interrupted  exposure  is  made.  This 
allows  measurement  of  leading-edge  to  leading-edge  distance  which 
limits  tail-off  errors  inherent  in  single  exposures.  The  measured 
particle  displacement  between  the  two  exposures  together  with  the 
interruption  time  (the  order  of  a  microsecond)  allows  calculation 
of  the  particle  velocity.  Because  of  the  size  of  the  zone  of 
interest  and  the  required  camera  location  to  keep  it  out  of  the 
hot  environment,  a  complicated  optics  system  is  required.  It  is 
also  anticipated  that  the  optics  will  include  the  capability  of 
photography  in  two  perpendicular  planes  and  thus  allow  a  three- 
dimensional  determination  of  particle  speed  and  trajectory. 

Particle  temperatures  both  before  and  after  impact  will  be 
measured  by  one  or  both  of  two  techniques.  One  technique  employs 
a  two-color  optical  pyrometer  which  views  the  stream  perpendicular 
to  the  flow  direction.  With  large-particle  loadings  for  which  the 
particles  appear  as  a  cloud,  the  pyrometer  will  read  out  the  cloud 
temperature  and  therefore  the  approximate  particle  temperature. 

For  light-particle  loadings,  its  utility  is  questionable  however. 
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The  possibility  of  making  a  tare  measurement  to  account  for  the 
pyrometer  response  to  the  stream  alone  and  of  correcting  for  the 
fraction  of  its  view  field  which  includes  particles  is  being  con¬ 
sidered.  Hope  is  small  here,  however.  The  other  method  is  appli¬ 
cable  to  light  particle  loadings  and  employs  a  technique  developed 
at  Vidya  (Ref.  26) .  A  double-image  photograph  of  the  particle¬ 
laden  gas  stream  and  of  a  calibrated  tungsten  filament  is  made. 
Temperature  is  then  determined  by  a  color  densitometer  analysis 
of  the  exposed  film  through  comparison  of  the  gas  stream  and  fila¬ 
ment  images.  The  use  of  this  technique  in  conjunction  with  the 
image  intensifier  camera  is  also  being  investigated. 

The  heat  flux  due  to  radiation  from  the  particles  will  be 
determined  by  radiometer.  In  some  cases,  the  heat  flux  into  the 
test  section  will  be  determined  from  the  measured  internal  tem¬ 
perature  distribution.  The  thermocouple  probes  and  data-reduction 
technique  will  be  similar  to  those  used  in  the  chemical  erosion 
program.  At  least  one  wall-temperature  history  will  be  obtained 
for  each  test  section. 

Particle  sampling  will  also  be  attempted.  A  number  of  tech¬ 
niques  are  under  study,  however,  none  seem  to  offer  a  reliable  and 
practical  answer.  The  problem  is  still  under  investigation. 

The  complete  definition  of  all  instrumentation  procedures  will 
be  completed  during  the  next  reporting  period  and  setup  begun. 

The  actual  start  of  testing  in  the  program  will  necessarily  depend 
on  the  time  required  to  get  all  instrumentation  set  up  and  operating. 
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TABLE  II 


CHEMICAL  EROSION  OF  GRAPHITE  IN  H,  C,  N,  0,  and  Cl  EXHAUST 
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•  1 05 1 E  +04 

2.89 

20.33 

1 • OOOE-OO 

.1178 

2 

3200.0 

10.0 

•  1466 

•8396E+03 

2.92 

20.03 

2. OOOE-OO 

.1178 

2 

3200.0 

25.0 

.2073 

.931 0E+03 

2.80 

20.34 

1 • OOOE-Ol 

.1178 

2 

3200.0 

25.0 

•  2073 

.931 0E  +  03 

2.80 

20.34 

1 . OOOE-OO 

.1178 

2 

3200.0 

25.0 

.1635 

•8152E+03 

2.61 

20*44 

2. OOOE-OO 

.1178 

2 

3200.0 

40.0 

.1986 

•8807E+03 

2.76 

20.3  1 

1 . OOOE-O! 

.1178 

2 

3200.0 

40.0 

•  1986 

•8807E+03 

2.76 

20.3  1 

1 • OOOE-OO 

.1178 

2 

3200.0 

40.0 

•  1  666 

•7965E+03 

2.76 

20.50 

2. OOOE-OO 

.1178 

2 

3200.0 

50.0 

.1950 

•8604E+03 

2.74 

20.30 

1 • OOOE-Ol 

.1178 

2 

3200.0 

50.0 

•  1950 

•8604E+C3 

2.74 

20.30 

1 .OOOE-OO 

.1178 

2 

3200.0 

50.0 

.1658 

•  7829E  +03 

2.74 

20.48 

2. OOOE-OO 

.1178 

2 

3600.0 

5.0 

.9195 

•2796E+04 

2.75 

26.18 

1 . OOOE-Ol 

.1178 

2 

3600.0 

5.0 

•  9195 

•2796E+04 

2.75 

26.18 

1 • OOOE-OO 

.1178 

2 

3600.0 

5.0 

•  9195 

.2796E+04 

2.75 

26.1  8 

2*  OOOE-OO 

.1178 

2 

3600.0 

1  0.0 

.7572 

•2494E+04 

2.73 

25.28 

1 .OOOE-Ol 

.1178 

2 

3600.0 

1  0.0 

.7572 

•2494E+04 

2.73 

25.28 

1 • OOCE-OO 

.1178 

2 

3600.0 

10.0 

.7572 

•2494E+04 

2.73 

25.28 

2*  OOOE-OO 

.1178 

2 

.0801 


Molecular 

Weight  20.36  19.36  19.54  11.90  11.96 


TABLE  IV 


CHEMICAL  EROSION  PROGRAM  TEST  CONDITIONS 


Test 

Point 

Mixture 

Material 

1 

1 

Graphite 

3500 

300 

(0.3) 

2 

1 

Silica  Phenolic 

3500 

300 

(0.3) 

3 

1 

Graphite  Phenolic 

3500 

300 

(0.3) 

4 

1 

Graphite 

2750 

300 

(0.3) 

5 

1 

Graphite 

2000 

300 

(0.3) 

6 

2 

Graphite 

3500 

300 

(0.3) 

7 

2 

Graphite 

3500 

300 

(0.3) 

8 

2 

Silica  Phenolic 

3500 

300 

(0.3) 

9 

2 

Graphite  Phenolic 

3500 

300 

(0.3) 

10 

2 

Graphite 

3500 

175 

(0.4) 

11 

2 

Graphite 

2750 

300 

(0.3) 

12 

2 

Silica  Phenolic 

2750 

300 

(0.3) 

13 

2 

Graphite  Phenolic 

2750 

300 

(0.3) 

14 

2 

Graphite 

2000 

300 

(0.3) 

15 

2 

Silica  Phenolic 

2000 

300 

(0.3) 

16 

2 

Graphite  Phenolic 

2000 

300 

(0.3) 

17 

3 

Graphite 

3500 

300 

(0.3) 

18 

3 

Graphite 

3500 

300 

(0.3) 

19 

3 

Silica  Phenolic 

3500 

300 

(0.3) 

20 

3 

Graphite  Phenolic 

3500 

300 

(0.3) 

21 

3 

Graphite 

3500 

175 

(0.4) 

22 

3 

Graphite 

2750 

300 

(0.3) 

23 

3 

Silica  Phenolic 

2750 

300 

(0.3) 

24 

3 

Graphite  Phenolic 

2750 

300 

(0.3) 

25 

3 

Graphite 

2000 

300 

(0.3) 

26 

3 

Silica  Phenolic 

2000 

300 

(0.3) 

27 

3 

Graphite  Phenolic 

2000 

300 

(0.3) 

28 

4 

Graphite 

3500 

300 

(0.3) 

TABLE  IV 


CONCLUDED 


Test 

Point 

Mixture 

Material 

Total 

Temperature 

Plenum 

Pressure 

(P£ 

3ia) 

29 

4 

Graphite 

3500 

300 

(0.3) 

30 

4 

Silica  Phenolic 

3500 

300 

(0.3) 

31 

4 

Graphite  Phenolic 

3500 

300 

(0.3) 

32 

4 

Graphite 

3500 

175 

(0.4) 

33 

4 

Graphite 

2750 

300 

(0.3) 

34 

4 

Silica  Phenolic 

2750 

300 

(0.3) 

35 

4 

Graphite  Phenolic 

2750 

300 

(0.3) 

36 

4 

Graphite 

2000 

300 

(0.3) 

37 

4 

Silica  Phenolic 

2000 

300 

(0.3) 

38 

4 

Graphite  Phenolic 

2000 

300 

(0.3) 

39 

5 

Graphite 

3500 

300 

(0.3) 

40 

5 

Graphite 

3500 

300 

(0.3) 

41 

5 

Silica  Phenolic 

3500 

300 

(0.3) 

42 

5 

Graphite  Phenolic 

3500 

300 

(0.3) 

43 

5 

Graphite 

3500 

175 

(0.4) 

44 

5 

Graphite 

2750 

300 

(0.3) 

45 

5 

Silica  Phenolic 

2750 

300 

(0.3) 

46 

r 

Graphite  Phenolic 

2750 

300 

(0.3) 

47 

5 

Graphite 

2000 

300 

(0.3) 

48 

5 

Silica  Phenolic 

2000 

300 

(0.3) 

49 

5 

Graphite  Phenolic 

2000 

300 

(0.3) 

50  ''l 

51 

52 

53  ( 

Conditions  as  found  desirable 

or  necessary  from  other 

54 

test  i 

-esults 

55 

56 

Test  nozzle  throat  diameter 


Figure 


Figure  2.-  Flow  diagram  of  wall-erosion  calculation. 


Figure  4.-  Measured  wall  temperature  axial  distribution,  JPL  Test  218 


Measured  wall  pressure  axial  distribution,  JPL  Test  218 


Figure  6.-  Calculated  energy  thickness  axial  distribution. 


Figure  7.-  Calculated  wall  heat  flux  axial  distribution. 


Figure  8.-  Comparison  of  results  with  experimental  data,  JPL,  Test  218 


Figure  9.-  Calculated  enthalpy  heat  transfer  coefficient  axial  distribution. 


Figure  10.-  Calculated  Stanton  number  axial  distribution 


Complete  Equilibrium 


Figure  12.-  Predicted  graphite  wall  erosion  in 
H,  C,  N,  O,  and  Cl  exhaust. 


Temperature,  K 

Figure  13.-  Pressure  and  mass  transfer  effects  on 
predicted  chemical  erosion  of  graphite  in  H,  C, 
N,  0,  and  Cl  exhaust. 


Figure  14.-  Surface  temperature  for  transient  heat  conduction 


Transient  conduction  in  the  hollow  cylinder,  back  wall  insulated 


[tYtTil 


Temperature ,  K 
(a)  Mixture  1 . 

Figure  16.-  Temperature  enthalpy  variation  for  test  gas  mixtures 


Temperature ,  K 
(b)  Mixture  2 . 


Figure  16.-  Continued 


Temperature,  K 
(d)  Mixture  4 


Figure  16.-  Continued 


Enthalpy,  Btu/lb 


Temperature 
Entropy ,  Bt 


Figure  17.-  Master  thermodynamic  charts  for  test  gas  mixtures 


Figure  17.-  Continued 


Figure  17.-  Concluded 
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Figure  18.-  Continued 
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(c)  Thermodynamic  properties  . 
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Figure  18.-  Continued. 


Hydrogen  }«— I  Oxygen  Wl  Mi^tur^Supply 


Figure  19.-  Gas  feed  and  metering  systems. 


Figure  20.-  Test  setup  for  arc-plasma  generator  calibration. 


Gas  enthalpy  determined  by  energy  balance, 

Btu/lb 
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Gas  enthalpy  determined  by  calorimeter 
(or  mass  discharge) ,  Btu/lb 

Figure  21.-  Comparison  of  gas  enthalpy  by  energy  balance, 
calorimeter  and  mass  discharge  methods  for  nitrogen  gas. 


Figure  23.-  Pressure-tap  locations  in  c  libration  nozzles 


Figure  24.-  Heat-transfer  calibration  nozzle 


Figure  26.-  Thermocouple  installation  in  test  nozzle 


Figure  28.-  Thermocouple  locations  and  analytic  temperature  distributions  in 

test  nozzles. 


Figure  28.-  Concluded 


Pressure  ratio 


(b)  0 .4-inch-diameter  throat 
Figure  29.-  Concluded. 


Figure  31.-  Surface  recession  as  determined  from  plenum  pressure 

Test  No.  795. 


Continued 


Figure  32.-  Nozzle  firing,  ARJ  graphite;  Test  No.  795. 


APPENDIX  A 


DISCUSSION  OF  A  STABILITY  CRITERION  FOR  THE 
FINITE  DIFFERENCE  FORMULATION  OF  THE 
HEAT  EQUATION 


It  is  stated  on  page  35  that  "the  general  condition  of 
stability  requires  that  the  coefficient  on  T  appearing  in 
Equation  (59a)  be  positive."  This  statement  requires  some  dis¬ 
cussion,  and  such  is  provided  in  this  appendix.  The  general 
condition  of  stability  to  which  reference  is  made  applies  to  the 
explicit  finite  difference  formulation  of  the  heat  equation. 

It  is  simply  stated  as  follows:  "The  present  temperature  of 
any  node  must  depend  on  the  past  temperature  of  that  node  in  a 
non-negative  sense."  Thus,  in  Equation  (59a),  the  coefficient 
of  T  must  be  positive  (or  zero) . 

The  authors  have  seen  no  mathematical  proof  of  this  "general 
condition,"  but  rather  arrived  at  it  by  a  heuristic  argument  made 
by  Schneider  (Ref.  18)  ,  and  have  developed  confidence  in  it  through 
successful  application  during  the  past  several  years. 

It  has  been  found  extremely  useful  for  determining  allowable 
time  increments  for  problems  complicated  by  such  things  as  geometry 
chemical  reactions,  and  convection  (at  the  surface  or  in  depth, 
as  through  a  char) .  In  the  simple  case  of  one-dimensional  con¬ 
duction,  this  requirement  leads  to  the  familiar  criterion  that 


a  Aa 
(Ax)2  ^ 


1 

2 


and  this  criterion  has  received  mathematical  proof  (see  Refs.  18 
and  28) . 

A  plausibility  argument  for  the  general  condition  of  stability 
which  the  authors  find  of  some  satisfaction  is  presented  below. 
Admittedly,  it  says  nothing  about  convergence  or  stability  of  the 
finite  difference  solution  per  se,  but  rather  is  based  on  the  cri¬ 
terion  that  the  formulation  satisfy  the  Second  Law  of  Thermodynamic 

i 

I 

1 


A-2 


and  thus,  in  this  sense,  be  "reasonable."  We  can  only  defer  to 
its  successful  use  to  relate  the  result  to  stability  of  the  solu¬ 
tion. 

Consider  a  simple  one-dimensional  body  with  the  temperature 
distribution  shown  in  the  sketch  below  at  time  Q. 


x 


The  customary  explicit  formulation  expression  for  the  temperature 
at  node  1  at  time  6  +  A6  is 

T  '  =  m  |T  +  T  +  (  —  -  2)  T  1  (A-l) 

l  12  3  m  l  I 


where 


ctAe 

(Ax)2 


For  simplicity  we  have  specified 


T 

2 


T  ,  so 

3 


T 

1 


m 


2T  + 
2 


(-  -  2)  T 
'  m  l  _ 


(A-2) 


Rearranging , 


=  T  -  2m  (T  -  T  ) 
1  12 


T  ’ 

l 


(A-3) 


Consider  now  the  three  possible  cases  for  m. 


(a)  (  £  -  2)  >  0 

In  this  case,  m  <  1/2,  and  from  Equation  (A-3)  and  exami¬ 
nation  of  the  sketch  we  have : 

T  1  <  T  and  T  '  >  T 
11  12 

(b)  (  I  _  2)  =  0 

m 

Here,  m  =  1/2,  and  T  1  <  T  ,  and  T  1  =  T 

11  1  2 

(c)  (  m  “  2)  <  0 

Here,  m  >  1/2,  and  T  '<  T  and  T  ’  <  T 

1^1  1^2 

Now  for  Ti ‘  <  T2 >  energy  must  be  transferred  from  node  1  to 
nodes  2  and  3  at  a  higher  temperature.  Since  this  violates  the 
Second  Law  of  Thermodynamics,  we  conclude  that  it  is  unreasonable 
to  specify  (1/m  -  2)  <  0.  And  in  fact,  in  practice,  such  specifi¬ 
cation  leads  to  large  oscillations  with  time  in  the  calculated 
temperature  of  node  1. 


ERRATA  TO  FIRST  QUARTERLY  PROGRESS  REPORT 


Page  6,  line  6 
Page  16, 

Page  31, 


Page  32, 

Page  35, 


Page  38, 
Page  39, 

Page  64, 
Page  66, 


The  number  0  should  be  replaced  by  the 
words  v30  percent  above".  .  . 

The  first  sentence  beginning  on  that  page 
should  read:  "in  theory,  any  value  of  T) 
could  be  used,  although  ..." 

The  last  two  equations  on  the  page  should 
be : 


k' 

2 


P  U  C», 
r  e  e  M 


k  1 

l 


k  p  M 
i^s  s 


^H20  peUeCM 


The  sentence  beginning  the  last  paragraph  on 
the  page  should  read:  "The  limiting  case  of 

k,'  =  k'  =  0  has  been  calculated." 

The  last  sentence  on  that  page  should  read: 
"Assuming  an  Arrhenius  extrapolation  of  the 
data,  that  is,  plotting  the  log  of  £he  coef¬ 
ficients  versus  1/T,  yields  at  3000  K 

mc  =  0.1  moles  carbon/min-cm2 

which  is  of  comparable  magnitude  to  observed 
rates. " 


The  second  sentence  on  the  page  should  read: 
"Patch  (Ref.  40)  has  studied  the  ..." 

Credit  for  the  derivation  of  the  equation  at 
the  top  of  the  page  should  be  given  to 
Dr.  R.  M.  Kendall  of  the  Vidya  Staff.  To 
the  author^1  knowledge  this  does  not  appear 
elsewhere  in  the  literature. 

The  first  sentence  below  Table  3.6  should 
indicate  the  molecular  weight  as  "11.92". 

The  R/M  term  in  Equation  (62)  should  read  as: 


Page  67 


Page  68, 


Page  70, 


The  tabulation  at  the  bottom  should  read: 
Simulated  Actual 


M 

Y 


1/2  inch 
18.60 
1.  17 


1  foot 
20 

1.  17 


The  equation  at  the  top  should  read: 
P„ 


sim 


=  0.510 


act 


The  second  line  after  that  expression 
should  indicate  "l53  psia". 


The  equation  should  read: 


sim 


act 


1. 075 


The  sentence  following  should  read: 

"For  this  case,  the  shear  stress  will  be 
some  7.5  percent  higher  in  the  laboratory 
simulation  than  in  the  actual  rocket  situ¬ 
ation.  " 

The  omitted  page  number  in  the  first  sentence 
of  'Summary'  is  page  55. 

Inadvertently  Table  I  was  omitted  from  the 
First  Quarterly  Progress  Report.  For 
completeness,  it  is  included  here. 
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